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Thoriated tungsten wires have been so heat-treated in 
gas as to develop variations in grain size, ranging in length 
from 44 to 3 cm. These variations produced marked 
changes in the thermionic emissive characteristics, such as 
rates of deactivation and of activation and duration of 
emissive life at normal operating temperature. The heat 
of diffusion remained unaffected. It was concluded that 
intergranular diffusion of thorium through tungsten is 


relatively too rapid to be measured and that the observed 
phenomenon depends purely on intragranular diffusion. 
This conception is in accord with the existence of a critical 
grain size for the realization of a maximum emissive life. 
Experiments are also described to demonstrate the 
spreading of thorium atoms over the surface of a single 
crystal of tungsten and to manifest the effect of slight 
strains in the lattice upon activation. 


XAMINATION of diverse lots of thoriated 
tungsten wires has shown that the values of 
the diffusion and activation constants may vary 
widely from the values cited by Langmuir.' These 
apparent discrepancies were finally found to be 
related to grain size and their explanation in that 
light has furnished fresh experimental evidence 
on the mechanism of diffusion of one metal 
within another, which is a development of the 
conception of Compton and Langmuir? that 
“diffusion is particularly rapid along crystal grain 
boundaries but must occur also from the interior 
of crystals to these boundaries and probably 
also through the surface layer of crystals.” 

The clue to this relationship, as well as the 
elimination of other possible causes, came from a 
study of the variation in emissive life, expressed 
as the length of time for the emission of a filament 
0.010 cm in diameter to drop to 20 percent of its 
initial value at a temperature of 2000°K. The 
construction of the kenotron and the general 
technique followed in all tests were the same as 


those described by Dushman, Rowe, Ewald and 


1]. Langmuir, Phys. Rev. 22, 357 (1923). 


°K. T. Compton and I. Langmuir, Rev. Mod. Phys. 2, 
138 (1930). 


Kidner.’ The life was found to be uniform for 
different filaments cut from the same wire lot. 
Seven different wire lots examined exhibited, 
however, a range in life varying from 140 to 1200 
hours. Chemical analysis showed variations in 
thoria content of the various lots ranging from 
0.9 percent to 2.8 percent but these were un- 
related to length of life. It was found, further- 
more, that the loss of thoria during life 
amounted to less than 10 percent for both the 
good and poor wires. The size of the thoria 
particles was found to vary somewhat but this 
also was unrelated to life. 

It then developed that the controlling factor 
lay in the heat treatment which the wire received, 
and that the effectiveness of this lay in the 
possibility of thereby increasing the grain size to 
a critical value. The life of & 200 hour wire, for 
instance, was increased seven-fold by a 1 min. 
flash at 3000°K in argon at 60 cm pressure. An 
incidental result of this is an increase in the 
concentration of reduced thorium by reason of 
the effect of a heavy gas in retarding the evapo- 
ration of a metal. The resistivity of thoriated 


3S. Dushman, H. N. Rowe, J. Ewald and C. A. Kidner, 
Phys. Rev. 25, 338 (1925). 
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wires which had been tested in 50 cm pressure of 
argon rose from the normal value for pure 
tungsten of 5.5 (10)~* ohms per cm’ to 5.8 (10)~*, 
denoting the liberation of thorium and its en- 
trance into solid solution with the tungsten. 
Such a change was absent in vacuum. This could 
have been only of minor importance upon life, 
however, for a wire made from a powdered 
mixture of tungsten with 1 percent of metallic 
thorium showed the abnormally short life of 90 
hours under normal treatment despite the fact 
that its temperature coefficient of resistance was 
20 percent below normal, demonstrating the 
retention of a large amount of metallic thorium 
in solid solution. The high temperature flash in 
argon had the same effect of increasing the grain 
size of this wire and this was accompanied by the 
same abnormal increase in life. Furthermore, the 
same favorable result came from a high tempera- 
ture treatment in hydrogen, too light a gas to 
retard appreciably the evaporation of thorium. 
In vacuum the treatment gave only. a slight 
improvement and in turn there was a greater 
resistance to grain growth. Even under the 
normal treatment a continuous improvement in 
life resulted during successive life tests, each 
preceded by the customary activation and ac- 
companied by a slight but continuous increase in 
grain size. 

Similar effects on grain size have been observed 
by Jeffries.‘ He has described the slight grain 
growth that occurs on long continued burning at 
low temperatures and has explained the effect of 
thoria in producing exaggerated grain growth 
when a wire is flashed on abruptly at tempera- 
tures above 2800°K. It was found possible by 
variations in this treatment to vary the grain size 
at will up to single crystals several cm long. 

The behavior of such wires with respect to 
emissive characteristics is described below. From 
these it developed that the fundamental phe- 
nomenon involved is one of intragranular diffu- 
sion. Its relation to grain size was found to 
explain the variations in life. 


DEACTIVATION 


An extended study was made of the deactiva- 
tion characteristics of different wires. The con- 


4 Zay Jeffries, Trans. Amer. Inst. Min. Eng. 60, 588 
(1919). 
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ventional method was followed—first flashing the 
filament for 3 minutes at 2800°K to form a 
supply of reduced thorium and then deactivating 
at the temperature in question. The course of the 
deactivation curve was found by interrupting 
this process at intervals and measuring the 
emission at 1450°K, a temperature so low that 
the status of the thorium with respect to diffusion 
and to evaporation was maintained unchanged. 

The diffusion coefficient, D, was calculated 
from the latter portion of the curve according to 
the second method employed by Langmuir,! 
involving the equation 


At 5.783 


in which @ is the fraction of the surface covered by 
thorium atoms at any time, ¢, in seconds and 8, is 
the limiting value of deactivation for a wire of 
radius Rin cm. 

The size of grains was determined from micro- 
scopic examination of the surface with the aid of 
the camera lucida at a magnification of 1000. 
For each filament the areas of about 400 grains 
were measured, the distribution was plotted and 
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Fic. 1. Effect of grain size on deactivation of thoriated fila- 
ments at 2400°K. A is fine grained, B is long grained. 
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DIFFUSION 


the effective grain radius taken as that of the 
area at the point of mid-distribution. 

An example is given in Fig. 1 to show the effect 
of grain size upon the deactivation curve at 
2400°K. Both filaments are from the same wire, 
(0.005 cm in diameter. A is the untreated sample 
whose effective grain radius is 5.4u, having the 
fine crystallization normally characteristic of 
thoriated wire. B had been given a special high 


TABLE I. Relation between grain size of thoriated wire and 
diffusion coefficient of thorium through tungsten at 2400°K. 


Grain Diffusion coeffi- Life at 
radius in » cient in cm? sec.~! 2000°K 
r Hours 
4.0 4.1 (10) 65 250 
4.8 2.5 160 
5.3 3.2 90 260 
6.2 1.5 58 3200 
6.3 2.2 87 3600 
6.9 1.3 62 18000 
7.3 1.3 69 5800 
8.2 0.9 61 9600 
8.6 1.0 74 18000 
1000 0.21 2800 
3000 0.01 50 


temperature treatment which served to convert 
it into a multiple grain wire comprising about 
seven grains across the area of the wire, each one 
measuring about 3 cm long. The value of D was 
found to be 1.7 (10)~*® cm? sec.—' for A, in the 
same range as the value reported for the same 
temperature by Langmuir, 3.6 (10)~*. For B, D 
was reduced about one hundred-fold, to 1.1 
(10)—", 

More comprehensive measurements were made 
on a lot of thoriated wire, 0.0102 cm in diameter 
and containing 1.5 percent thoria. Samples were 
so heat treated as to cover a wide range in grain 
size. In the last two cases the grains filled entirely 
the cross section of the wire and had the average 
length noted. The results are given in Table I. 

Table I shows a distinct trend of D toward 
lower values as the radius of the grain increases. 
The significance of this is that it is not the 
radius of the wire but the radius of the grain 
which should be incorporated in the formula for 
the diffusion. The original derivation presumed a 
certain concentration gradient throughout the 
area of the wire as.a whole as expressed in the 
equation, G= (2C,,/R)e—*-78"/*", The present re- 
sults denote that it is rather the area of the grain 
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which is the important entity and that the 
concentration gradient should be with reference 
to the distribution of thorium within each grain, 
with the term R? referring to the radius of the 
latter. This is brought out still more clearly in the 
3d column of Table I, where a revised version of 
the diffusion coefficient is expressed as the 
product of D, the conventional coefficient, and 
r’, the square of the effective grain radius 
characterizing each filament. To within the limits 
of error in measuring grain size, the product 
remains constant. Such a relationship denotes 
that the rate of diffusion as measured is that 
characteristic of diffusion within a grain and that 
the rate of diffusion along grain boundaries is 
relatively too rapid for measurement. 

Saul Dushman and L. R. Koller have studied 
in this laboratory the deactivation of thoriated 
wires for differences in diameter. Their un- 
published results are also in accord with the view 
that grain size is a predominant factor in de- 
termining the magnitude of the diffusion con- 
stant. 

Additional support is given to this view by 
measurements on the temperature coefficient of 
diffusion. In Table II are given values of D for: 
several different temperatures as determined for 
filaments of different grain size. For comparison, 
the Langmuir values are given in the lowest 
horizontal row. 


TABLE II. Diffusion coefficient of thorium through tungsten 
at various temperatures. 


Grain Diffusion coefficient in cm? sec.~! 
radius in yu D 
r 2400° 2300° 2200° 2000° 
5.3 3.19(10)-® 1.48(10)~ 5.8 
7.3 1.34(10)-* 6.6 (10)-” 2.18(10)-" 
3000 1.13(10)-" 2.26(10)-"3 


Langmuir 3.02(10)~* 1.23(10)-" 4.90(10)- 


The logarithms of these values of D are plotted 
against the reciprocals of the temperature in Fig. 
2. The curves are displaced relatively to one 
another because of the differences in grain size 
involved. The slope of all three, however, is 
identical and the same in fact as that found by 
Langmuir, corresponding to a heat of diffusion of 
94,000 calories per g-atom. If both types of 
diffusion took place at measurable rates, then 
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Fic. 2. Diffusion of thorium as function of temperature 
and grain size. A, 5.34 radius. B, 7.34 radius. C, 3000u long. 


corresponding differences in the heat of diffusion 
would appear in experiments covering such a 
wide range in grain size and one or the other 
type of diffusion would predominate at the 
extremes. 

In the last column of Table I are given values 
for the emissive life at a temperature of 2000°K 
to bring out the effect of grain size upon it. The 
life evidently begins to become appreciable only 
at a grain size whose radius is about 6u. Then it 
increases somewhat at least up to a grain size 
whose radius is 9u. Thereafter it decreases, 
eventually becoming negligible for sufficiently 
long grains. 

This is explained from the assumed relation 
between intergranular and intragranular diffusion 
and from their relation in turn to the phe- 
nomenon of induced evaporation. Whatever the 
mechanism of the latter may be, it seems to be 
unquestioned that a layer of thorium more than 
one atom thick is relatively unstable and that the 
outer layers will evaporate more readily than 
those atoms which are bound directly to tung- 
sten. After the high temperature flash associated 


with activation, a filament has within it a certain 
concentration of reduced thorium. This is being 
continually added to, although only in small 
amounts, by the thorium that is reduced at the 
normal operating temperature. It is upon the 
initial concentration that dependence must be 
placed primarily to maintain such a gradient 
that the surface of the wire may be kept covered 
with thorium and a decrease in emission avoided. 
The maximum life would result when the loss by 
normal evaporation is just compensated for by 
diffusion of an equal amount of thorium from 
within the filament. The rate of this diffusion is 
given by the product DG, the diffusion constant 
and the concentration gradient. Only within 
certain critical values of D will there be an 
approach to the ideal compensation desired. If 
D is too small, which is the condition in a long 
grain filament, then the diffusion is obviously 
insufficient to maintain the full thorium coating 
on the surface of the wire. If the grains, on the 
other hand, are so fine as to give rise to too large 
a value of D, an equally undesirable condition 
exists, for the total amount of thorium diffusing 
is then in excess of the normal loss by evaporation 
and this excess evaporates from the surface of the 
wire so rapidly, through the phenomenon called 
induced evaporation, that it is virtually wasted. 
Such a loss lowers in turn the concentration 
gradient of reduced thorium which thereby falls 
relatively soon to such a low value that it fails to 
maintain an amount of diffusion that would be 
the equivalent of the amount normally evapo- 
rating. 


ACTIVATION AND SURFACE SPREADING 


The rate of activation, k, is given by the 
expression d0/dt=k(1—6). At a temperature of 
2050°K the evaporation of an adsorbed film of 
thorium is so extremely slow that k is propor- 
tional to the product DG (Langmuir'). One 
should expect therefore that the differences that 
have been found in the diffusion constant, as 
defined by D, should be accompanied by similar 
changes in k. Such differences in k have already 
been noted in unpublished work carried out in 
this laboratory by D. M. Dennison. 

This relationship is manifested by the present 
experimental results, as is evident from the 
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definite trend shown in Table III. The determi- 
nation of k was in each case preceded by a 3 
min. flash at 2800°K to bring G to a constant 
value. 

These results suggested a means of determining 
whether thorium atoms can spread over the 
surface of a filament. In the light of the present 
conception, thorium is diffusing to the surface 
throughout the length of the filament. The 
diffusion outward through a grain boundary will 
convey, however, many more thorium atoms to 
the surface at the intersection of two grains and 
therefore give rise to a higher concentration along 


TaBLe III. Relation between grain size and activation 
constant at 2050°K. 


Grain radius Activation constant 
k 


4.0 10.9 (10)-* 
4.8 17.2 
10.0 
6.3 8.9 
6.9 4.2 
5.9 
8.2 4.6 
8.6 3.3 
3000 3.2 


the boundary than can take place on the surface 
of a grain as the result of direct diffusion only 
through the grain itself. This is by reason of the 
large internal volume of filament that feeds 
into a grain boundary. It was aimed to find if the 
high thorium concentration formed at grain 
boundaries led to a spreading of thorium over the 
surface of adjacent grains. The experiment was 
carried out by making up a filament in the form 
of a U, each leg being 7 cm long joined by a 1 cm 
length at the bottom. The legs were converted 
into grains about 2 cm long by flashing on 
separately at 3000°K in hydrogen and raising 
the temperature to 3300°K in the course of one 
minute. The 1 cm length at the bottom was 
rendered finely granular with an average grain 
radius of 4u by raising in temperature from 1200° 
to 3000°K in the course of 5 min. When this 
filament was made into a kenotron, its rate of 
activation at 2050°K was found to be 5.3 (10)-*. 
The emission as read at 1500°K rose to a maxi- 
mum value of 1780ua as compared to a value of 
5000ua customarily found for a filament of the 
same total length, 15 cm, but entirely fine 


grained. This corresponds to a value of 6 for the 
entire filament of 0.88. The fine grain portion of 
the filament represented 6.7 percent of the total 
area and if the covering of the surface by 
thorium had been confined to it, the emission 
would have risen only to 334ua. The validity of 
the experiment was established by cutting off 
the fine grain portion of the filament and re- 
mounting the remainder in a length of 11 cm. 
The activation constant was thereupon found to 
be reduced to 1.6 (10)-5. At the end of 5 hours, 
the time required to reach maximum emission 
when the fine grain portion was attached, its 
emission had risen only to 32ua. The results give 
a clear demonstration, therefore, of the depend- 
ence of the activation constant upon grain size 
and also of the occurrence of spreading of 
thorium atoms over the surface of tungsten. 

Gehrts® has derived the Langmuir equation for 
activation on different assumptions—namely, 
that activation depends primarily upon the 
speed with which thor'um atoms spread over the 
surface of a crystal rom grain boundaries as 
sources and that induced evaporation is absent. 
In certain special cases, as in that just described 
of a single crystal adjoining a mass of fine crys- 
tals, the speed of activation of the composite 
filament depends, of course, upon the speed of 
spreading and it would be possible to perform the 
experiment in such a way that the speed of 
spreading over the single crystal could be meas- 
ured separately. For filaments of uniform grain 
character, however, the observations given above 
on the dependence of deactivation and life upon 
grain size would seem to involve diffusion as the 
fundamental phenomenon. 


EMISSION OF SINGLE CRYSTALS AND STRAIN 


Clausing® concluded that intragranular diffu- 
sion is altogether negligible because of his 
inability to develop any activation in Pintsch 
single-crystal filaments. We found, likewise, that 
Pintsch filaments could not be activated to any 
appreciable emission at the normal activating 
temperature of 2050°K or higher. In the light of 
the present experiments this does not mean that 
intragranular diffusion is negligible but that it 


5 Gehrts, Zeits. f. tech. Physik 12, 66 (1931). 
® Clausing, Physika 7, 193 (1927). 
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can advantageously be measured at this tempera- 
ture only on a multigranular wire. At 2050° the 
diffusion in a single crystal brings so little 
thorium to the surface in unit of time that the 
most of it is lost by normal evaporation and no 
opportunity given for accumulation of any per- 
ceptible amount. We were able, however, to 
activate the Pintsch single-crystal wire to its 
full emission by heating for 600 hours at 1750°K, 
a temperature so low that evaporation does not 
enter as a disturbing factor. 

From tests of Pintsch filaments in incandescent 
lamps there was evidence that mechanical strains 
could be introduced that were insufficient to 
cause recrystallization but which were neverthe- 
less severe enough to weaken the lattice structure. 
Filaments 0.004 cm in diameter were coiled on a 
0.025 cm mandrel and then uncoiled, so that 
their final form was once more that of a straight 
filament. When these were tested at 2600°K in 
vacuum lamps, their life was found to be reduced 
to 20 percent for the reason that they failed by 
burn-out after losing only 9.6 percent in weight. 
Untreated filaments continued burning until they 
had lost 23.6 percent in weight. Similar tests on 
drawn wire failed to show any differences due to 
coiling. The Pintsch filaments which had been 
coiled retained their single-crystal structure so 
far as outward appearances showed. Microscopic 
examination at 1000 magnification after burning 
showed the uniform surface of a single crystal, 
free from any grain boundaries. Burning the 
filaments in argon developed a hexagonal cross 
section, which was uniform throughout their 
length. This is a very delicate test, as the 
presence of a separate grain shows up distinctly 
by the abrupt change in orientation of the flat 
crystal faces. The conclusion therefore was that 
the coiling had not led to recrystallization but 
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only to the shattering of certain transverse 
lattice sections and that premature burn-out had 
occurred at one of these points. 

In order to find to what extent such effects 
might simulate grain boundaries, a_ similar 
Pintsch filament was given the same coiling 
treatment and tested in a kenotron. It did in 
fact behave like a fine grain filament activating 
at the normal rate of 9.7 (10)~* to full, normal 
emission. Despite this behavior which was typical 
of that of a fine grain filament, microscopic 
examination of the surface failed here likewise to 
reveal any departure from a single crystal, 
indicating that whatever sections were serving to 
simulate grain boundaries in furthering the 
diffusion of thorium were of atomic dimensions. 
They can perhaps be regarded as slight displace- 
ments of the lattice, insufficient to introduce any 
marked change in orientation. The effect is of 
similar nature to certain colorations observed by 
Przibram’ in single crystals of sodium chloride 
under radiation from 8- and y-rays from radium. 
When one portion of the crystal was very 
slightly deformed by pressure from a cube face, 
diagonal planes of coloration were produced 
throughout the underlying portion at an exposure 
insufficient to produce a general coloration. 
These planes coincided with the direction that 
would be expected as the result of a mechanical 
displacement. 

After this filament had been burned for 550 
hours at 2000°K, it was reactivated at 2050°K. 
The value of K was found to have fallen to 
1.4 (10)~, indicating a partial healing of the 
hypothetical displaced section of the lattice. 

In closing we wish to express our appreciation 
of the continued interest and encouragement 
shown us by Dr. Saul Dushman. 

7 Przibram, Zeits. f. Elektrochemie 38, 490 (1932). 
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Experiments on the Contact Potential of Zinc Crystals 


W. A. ZisMAN AND H. G. YAmins, Research Laboratory of Physics, Harvard University 
(Received November 7, 1932) 


T is well known that the experimental study of 
the Volta effect of metals in gases is difficult. 
Two facts have been established. Firstly, differ- 
ent methods of preparing the metal surfaces such 
as the use of various types of abrasives or cutting 
tools caused the contact potential difference to 
vary as much as 100 percent. Secondly, the po- 
tential difference varied with the time after prep- 
aration of the surface, apparently due to gas ad- 
sorption and corrosion of the metal. 

Despite the above difficulties every few years 
experiments appear in the journals repeating the 
same sort of experiment and leaving everyone 
more firmly convinced than ever of the futility of 
studying the Volta effect except under high 
vacuum conditions. Of course, research on thor- 
oughly degassed metals in high vacuum is neces- 
sary to obtain intrinsic contact potential differ- 
ences. Yet further work on the Volta effect when 
adsorbed gases are present in large amounts is 
interesting because of its intimate connection 
with the study of adsorption and corrosion. In 
this connection the vibrating condenser method! 
used by us is ideal because of the ease of taking 
measurements in rapid succession after prepara- 
tion of the surface. 

It has never been ascertained what share of the 
variability with methods of cutting and polishing 
is due to the alterations in the structure of the 
surface of the metal. The cleaved surface of a 
single crystal is the nearest to a reproducible one 
that we can obtain, and hence gas should be ad- 
sorbed on such a surface in the most repeatable 
manner possible. Consequently, the contact po- 
tential difference between a single crystal cleaved 
in air and some standard metal should be re- 
producible. 

Single crystals of pure zinc obtained from the 
New Jersey Zinc Company were cleaved along 
the basal plane, and the contact potential differ- 
ence with a standard electrode was observed as it 
varied with time. The standard electrode was a 
gold plate well aged in air. Measurements of the 


potential differences were made with the vibrat- 
ing condenser method of Zisman.' In our appa- 
ratus the vibrator plate was fastened to the dia- 
phragm of a telephone receiver driven by an 
audiofrequency vacuum tube oscillator. 

In Fig. 1 is given a typical graph of the change 
of the potential difference with time after cleav- 
age of the crystal. Some ten such cleavages were 
made, and the time curves were found to be re- 
producible to 1 percent. This is a vast improve- 
ment over any other work known to us on the 
Volta effect of zinc in air. For instance, if the 
surface of polycrystalline zinc is prepared by 
turning in a lathe with a clean tool, the time 
curves obtained for successive turnings vary by 
10 percent. If the surface is scraped, the variation 
is as much as 30 percent. 

A comparison of the time curves for turned 
polycrystalline zinc with those for a cleaved 
crystal reveal that the potential differences are 
about 40 percent lower for the former. For ex- 
ample, a number of runs on polycrystalline zinc 
netted a mean potential difference with gold of 
0.85 volts five minutes after turning, whereas the 
potential difference of a cleaved crystal at the 
same time was 1.22 volts. On the other hand, 
when the basal cleavage plane of a single crystal 
was turned in a lathe, the e.m.f.-time curves were 
precisely like those for the turned surface of 
polycrystalline zinc. This suggested that the 
lower values for the turned surfaces as compared 
with cleaved surfaces might be associated with 
the relatively greater roughness of the former. It 
is well known that roughness of surface is be- 
lieved to bring condensation due to the lower 
vapor pressure of a liquid in a small cavity as 
compared with that at a flat surface. However, 
when the single crystal was scraped in a lathe 
with a razor blade, the time curves showed that 
the potential difference was 0.65 volts five min- 
utes after scraping, which is even lower than for 
the turned surface, though the scraped surface 


''W. Zisman, Rev. Sci. Inst. 7, 367 (1932). 
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Fic. 1. Contact potential difference between a zinc crystal and a standard gold plate. 


was smoother than the turned. Thus, the lower 
value of contact potential difference cannot be 
associated with greater roughness. 

The turning or scraping of a metal like zinc 
involves a rolling and also a tearing action. It is 
well known? that rolling zinc causes about 70 per- 
cent of the small crystals in the surface to orient 
with their hexagonal axes making an angle with 
the metal surface of less than 45°, i.e., rolling 
causes a departure from the orientation existing 
in the metal before the surface was cut or 
scraped. Now, if the small crystals of zinc have a 
greater contact potential difference with gold on 
the basal plane than on the prism faces, then pre- 
cisely the observed effects would be predicted. 
Turning a single crystal would cause a drop in the 
contact potential difference because it caused the 
cleavage surface to be broken into small crystals 
no longer oriented in the original direction with 
the hexagonal axes perpendicular to the surface. 
Scraping would make the potential difference 
drop even more, because it involves more of a 
rolling effect than does turning. The same order 
of magnitude of potential difference would be 
expected whether we turned a single crystal, a 
polycrystalline aggregate composed of very small 
crystals, or one containing very large crystals. 


2 Wien-Harms, Handb. d. Physik 5, 79 (1932). 


Such was found to be the case. Successive turn- 
ings using the same type of tool and tool adjust- 
ment do not produce precisely the same distribu- 
tion of orientations among the crystals each time, 
and hence the potential difference for turned 
surfaces would be expected to vary from run to 
run, as has been observed. 

In accordance with these ideas, one would not 
expect any such variability with successive turn- 
ings when a metal like polycrystalline copper is 
used instead of zinc, for the copper crystals are 
each cubic, and therefore isotropic, whereas zinc 
is hexagonal and anisotropic. The anisotropy of 
zinc with respect to electrical resistivity and 
thermoelectric effect is well known. Thus, it is not 
surprising that it is also anisotropic with respect 
to the Volta effect. 

Successive turnings of polycrystalline copper 


_ina lathe produced e.m.f.-time curves which were 


identical to within 34 percent. Furthermore, un- 
like zinc, the scraped copper surface behaved 
precisely like the turned one. These two results on 
copper, reproducibility with successive turnings 
and identity of results of turning and scraping 
agree very well with what would be predicted by 
the above idea of anisotropy of the Volta effect. 

In connection with this idea of anisotropy re- 
cent experiments of Nitzsche’ are of interest. He 


3 A. Nitzsche, Ann. d. Physik [5], 14, 463 (1932). 
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found from a study of the photoelectric work 
function of cleaved surfaces of zinc crystals that 
the work function for the basal cleavage plane 
was 0.2 volts greater than that for the secondary 
cleavages planes perpendicular to the basa! plane. 
We attempted to obtain zinc crystals that could 
be cleaved in these secondary planes as well as 
the basal plane, but so far have been unable to do 
so. Further experiments in this direction are 
planned, for it would allow us to obtain more 
definite evidence of an anisotropy of the Volta 
effect in zinc. 

Some experiments on the use of abrasives like 
emery or of polishing agents like rouge in smooth- 
ing the surfaces of metals for the study of contact 
potential differences may be of interest here. 
Experiments dating as far back as 1880 demon- 
strated that the contact potential difference va- 
ried for surfaces rubbed with various types of 
abrasives like emery paper and sandpaper. As in 
the above case of a turned surface, the evident 
interpretations made then and continued until 
quite recently have been that according as one 
roughened the surfaces of the metal differently, 
so the adsorbed gas layer was altered, and hence 
the potential difference changed. In other words 
the variability obtained with abrasives was as- 
cribed to the difference they caused in the rough- 
ness of the metal. In 1929 Ende‘ reported experi- 
ments on the Volta effect in air using two pieces 
of metal from the same bar such as brass or nickel. 
The two like metals were polished with emery or 
else scraped or turned in a lathe. The first treat- 
ment produced contact potential differences be- 
tween the two like metals as great as 0.5 volts, 
whereas in the last two the resulting potential 
differences were always less than 0.05 volts. He 
suggested that emery left imbedded in the metals 
gave rise to a false potential difference of as much 
as 0.5 volts. 

We found that polycrystalline zinc or a cleaved 
crystal when polished with emery gave rise to a 
potential difference with gold of about 1.30 volts. 
When the surface was washed with water and 
dried the potential difference dropped to 1.10 
volts. This decrease was not due to any effects of 
the water, for when a polycrystalline surface free 


from emery was well washed with water, the 


‘ W. Ende, Phys. Zeits. 30, 477 (1929). 


potential difference after drying was within 0.01 
volts of the value before washing. Thus, the 
polishing of turned polycrystalline zinc with 
emery paper raised the potential difference from 
0.85 volts to 1.30 volts, because the polishing im- 
bedded emery in the zinc. Washing with water 
lowered the potential difference to 1.10 volts, be- 
cause some but not all of the emery was removed. 
In the same manner, polishing the zinc surface 
with rouge caused the potential difference to drop 
to 0.6-0.7 volts. After washing with water it 
rose to 0.8-0.9 volts, due to the removal of some 
of the rouge. These results therefore demonstrate 
the correctness of Ende’s suggestion. A number 
of experiments have been published in the past 
decade in which emery paper was used on metals 
which were later introduced into high vacuum 
systems for measurements of contact potential 
differences or thermionic work functions. The 
above experiments show the use of emery to be an 
unexpected source of complications. 


CONCLUSION 


Experiments were made proving for the first 
time that it is possible to obtain a reproducible 
contact potential difference on zinc in air at 
ordinary pressures, provided the cleaved surface 
of a single crystal is used. It is suggested that the 
lack of reproducibility of the Volta effect of 
turned or scraped surfaces of zinc is due to an 
anisotropy in the zinc crystals with respect to the 
Volta effect, and not due to condensation phe- 
nomena accompanying the roughness of the sur- 
face. Confirmatory evidence is found in the 
reproducibility of the contact potential difference 
of turned and scraped surfaces of copper. Ende’s 
suggestion that emery imbedded in the metals 
gives rise to false potential differences is demon- 
strated to be correct. It is found that emery in- 
creases the potential difference between zinc and 
gold by about 0.4 volts, and rouge decreases it by 
about 0.5 volts. 


We wish to thank the New Jersey Zinc Com- 
pany for furnishing us with the two single crystals 
of pure zinc used for these experiments, and to 
express our appreciation to the department of 
physics of Harvard University for their courtesy 
in allowing us to perform these experiments in the 
new research laboratory. 
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On the Piezoelectric Properties of Tourmaline 


GERALD W. Fox AND Morris UNDERWOOD, Physics Laboratory, Iowa State College 
(Received October 13, 1932) 


Oscillating piezoelectric plates of different frequencies 
have been made from California and South African 
Tourmaline. These plates have a frequency response of 
approximately 3770 ke per millimeter of thickness. The 
temperature coefficient for both specimens was found to be 
negative and is about 35.5 parts per million per degree 


centigrade for the African specimen and 38.1 parts per 
million per degree centigrade for the California crystal. 
Comparison of the tourmaline plates with a Y-cut quartz 
plate of the same size and frequency shows the tourmaline 
oscillators decidedly inferior on the basis of power output. 


INTRODUCTION AND STATEMENT OF THE PROBLEM 


INCE Cady’s! original discovery, the use of 

quartz plates to stabilize the frequency of 
vacuum tube oscillators has become a universal 
practice. Intense experimental investigation? has 
led to applications of great practical importance 
as well as to discoveries that yield considerable 
information concerning the general nature of the 
crystal state. According to F. Pockels,* twenty of 
the thirty-two possible crystal classes should 
show piezoelectric properties. However, simple 
exhibition of piezoelectric properties does not 
mean that a crystal will be useful for frequency 
control. This is well evidenced in the case of 
Rochelle salt which is highly piezoelectric but is 
mechanically so weak that it is entirely unsuited 
for use as a resonance oscillator. 

In the entire list of known piezoelectric crys- 
tals, there are few of value for frequency control. 
Aside from quartz and tourmaline and possibly 
boracite, all appear worthless for one reason or 
another. 

The present investigation of the piezoelectric 
properties of tourmaline proposed to determine 
its frequency constant and temperature coeffi- 
cient, and compare it so far as possible with 
quartz. 


DESCRIPTION OF THE CRYSTAL 


Crystallographically, tourmaline belongs to the 
trigonal system. It is an example of class 20 


!W. G. Cady, Proc. Inst. Rad. Eng. 10, 83-114 (1922). 

2 P. Vigereaux, Quartz Oscillators and Resonators. Excel- 
lent Bibliography, H. M. Stationary Office, London. 

3 F. Pockels, Winkelman’s Handbuch der Physik, Vol. 2. 
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(ditrigonal pyramidal) and is of the type ditrig- 
onal polar, having one trigonal axis and three 
symmetry planes intersecting in it. Its chemical 
constitution is somewhat variable though its 
formula is sometimes written as (HeNasF4BeAls- 
SijOx¢3). A diagram of the crystal model is shown 
in Fig. 1. The Z axis is the polar axis as well as the 
optic axis. Pressure applied along this axis results 
in a separation of charges at the ends. 


+ 


TT | 


Fic. 1. Crystallographic model of tourmaline. 


Tourmaline is a fairly common mineral. Gem 
quality crystals found in the United States come 
principally from Oxford County, Maine, and San 
Diego County, California. Most of the world’s 
supply, however, comes from Brazil and South 
Africa. Crystals used as piezoelectric oscillators 
must be especially perfect and even those which 
pass for gem quality are usually not good enough. 
Minute flaws and cracks throughout the body of 
the crystal are very common and, although such 
imperfections are tolerated in other uses, for 
oscillating plates, the section must be entirely 
free from them. The crystals occur in a great 
variety of colors, often with different colors ap- 
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pearing in the same specimen. Tourmaline is 
harder and heavier than quartz, hardness being 
7.5 and specific gravity about 3. On the gem 
market, prices for the raw crystals vary some- 
what, but average around one dollar per carat for 
the best quality. Unfortunately, the cross- 
sectional areas of these crystals are small; one 
having a cross section of one square centimeter 
would be considered large. 


PREPARATION OF PLATES 


From a large number of specimens, two were 
chosen as having possibilities for use as oscillators. 
One of these was from California and was ap- 
proximately three cm long and about eight mm 
in largest diameter. It shaded from pink at one 
end to light green at the other. Except for a six 
mm section near the center, this crystal was 
flawed for its whole-length. The other, a South 
African specimen, was about the same size. It was 
deep green in color and was only slightly flawed 
near one end. With a lapidary’s saw, fifteen sec- 
tions perpendicular ‘to the Z axis (Fig. 1) were 
cut from these two crystals. These pieces varied 
from 3 mm to 0.5 mm in thickness. Finishing the 
plates is a tedious process, at least when done by 
hand. It consists of grinding down the rough 
blanks by rubbing them on a plate glass slab 
with water and fine carborundum powder. The 
blanks are so small that it is difficult to hold them. 
The entire process is similar to the preparation of 
quartz plates except that it seems to require more 
care. In the case of quartz, it is well to have the 
faces parallel to within 0.01 mm, but this accu- 
racy is not good enough for tourmaline. Likewise, 
the edges of the plates need to be rounded off very 
carefully. The most minute unevenness or chip 
appears enough to prevent oscillation. Because of 
the small cross sections, no attempt was made to 
cut round disks from the blanks. Their finished 
shape corresponded to that of the crystal cross 
section except that the corners were rounded off. 
Fig. 2 shows two of the finished oscillators. The 
dark one is from the African crystal; the light 
one, from the California specimen. 


WAVE-LENGTH AND FREQUENCY CONSTANT 


The curve of Fig. 3 shows the wave-length to 
which the plates were resonant plotted against 


their thickness in mm. The points fall closely on a 
straight line, each point representing a separate 
finished plate. Both kinds of tourmaline used are 
represented, so that one would conclude the color 


No. 


Fic. 2. Finished tourmaline oscillators. No. 3 from S. 
African crystal, No. 10 from California specimen. 


is not an important factor. The equation of the 
curve is very closely: \=79.5 m where X is the 
wave-length and m is in meters per millimeter. In 
terms of frequency: 2= 3770 kilocycles per milli- 
meter. Recently Straubel* announced the use of 
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Fic. 3. Wave-length thickness curve for tourmaline. Plates 
from both specimens are included in the curve. 


tourmaline plates for ultra high-frequency con- 
trol. His plates were made for him by the firm of 
Carl Zeiss and from his article one would con- 
clude that it was quite simple to make tourma- 


4H. Straubel, Phys. Zeits. 32, 937-941 (1931). 
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line plates having a fundamental wave-length of 
five meters, i.e., fundamental frequency of 56 
megacycles and even higher. At the present time 
in this country, certain firms have such plates on 
the market and in a recent article, their use was 
described. From the curve (Fig. 3), one can 
readily estimate the thickness of a five-meter 
plate. A twenty-meter tourmaline plate is only 
slightly thicker than a quartz plate of the same 
frequency. Both are so thin as to be very weak 
mechanically. A five-meter tourmaline plate 
would be thinner than the paper on which this is 
printed and, as such, would be a very delicate af- 
fair indeed. From our experiments, it seems 
doubtful if plates having a fundamental of this 
high value can be successful practically. 


TEMPERATURE COEFFICIENT 


The temperature coefficient was determined 
for both specimens of tourmaline. Two plates cut 
from each crystal were carefully ground until 
their frequencies were within about two hundred 
cycles of each other. All four plates had frequen- 
cies close to 3500 kc. The conventional oscillator 
circuit was used in which the crystal is connected 
between the grid and ground, in parallel with a 
resistance for obtaining the proper bias. Two 
separate units with pentode tubes (UX 247) pro- 
vided plenty of output even though the plate 
voltages were kept low (100 volts) to avoid un- 
necessary heating. Only a few minutes were re- 
quired for steady conditions, and the beat note, 
heard on a high-frequency receiver with loud 
speaker, would remain constant for long periods. 
The room temperature was found not to vary 
appreciably over the one and a half hour period 
necessary for a run. A heavy brass plate with a 
light upper electrode provided the mounting for 
the constant temperature plate. The other 
mounting was similar except that the whole as- 
sembly was encased in a Pyrex tube and sur- 
rounded by a heater coil. A thermometer reading 
to 0.1°C was hung along side this tube for tem- 
perature determinations. 

The beat frequency was measured by a siren 
disk driven by a smoothly variable d.c. motor 
with heavy duty storage battery supply. A 
Weston voltmeter tachometer measured the mo- 


5 James Millen, Radio News 14, 206-208 (1932). 
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TABLE I. Frequency variation with temperature for plates 3 
and 5 (African tourmaline). 


Temperature Siren speed Frequenc 
(0°C) r.p.m. No. holes (calc.) J 
27.5 645 54 580 
29.4 410 54 369 
35.2 576 54 467 
38.8 1290 54 1161 
41.2 995 96 1592 
46.2 1275 96 2040 
49.0 1560 96 2496 
41.3 995 96 1592 
37.1 1115 54 1004 
35.0 720 54 648 
27.6 450 48 350 


TABLE II. Frequency variation with temperature of plates 2 
and 10 (California tourmaline). 


Temperature Siren speed Frequency 
(0°C) r.p.m. No. holes (calc.) 
26.4 795 96 1272 
30.5 595 48 476 
34.6 0 0 0 
39.2 905 48 724 
40.8 1425 48 1140 
44.1 1075 96 1720 
49.1 1225 96 2464 
45.7 540 96 1960 
38.6 0 96 428 
33.0 615 0 0 

30.1 1215 48 


492 


tor speed directly. This scheme proved entirely 
satisfactory for the precision expected and it was 
quite easy to determine zero beat between the 
siren note and the beat frequency of the crystai 
oscillators. Tables I and II give typical sets of 
data starting at room temperature and going to 
50°C and back. The curves of Fig. 4 are these 
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Fic. 4. Temperature frequency curves for tourmaline. 
The upper curve represents the African crystal, the lower 
one the California oscillators. 
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data plotted. The frequency variation is evi- 
dently linear over the range investigated and, as 
calculated from the upper curve (African tourma- 
line), the temperature coefficient is approximately 
35.5 parts per million per degree and is negative. 
For California tourmaline (lower curve), the co- 
efficient is also negative and, as calculated is 38.1 
parts per million per degree. 


COMPARISON OF TOURMALINE WITH QUARTZ 


In order to compare the output of quartz and 
tourmaline, it was necessary to make a quartz 
plate of about the same dimensions as the tour- 
maline plates. Since in practice, the quartz disks 
are nearly always larger than a centimeter in 
diameter, this meant a disk of about half the 
customary size. A Y-cut disk of 7 mm diameter 
was made having a fundamental frequency of 
3580 ke. With one of the pentode oscillator cir- 
cuits mentioned previously, observations were 
made on the power dissipated in an auxiliary 
coupled circuit consisting of a high-frequency am- 
meter in series with a loop made of nichrome rib- 
bon. The auxilliary circuit was coupled closely 
enough to the driver circuit to load it as much as 
possible without throwing it out of oscillation. 
The coupling was kept the same throughout the 
experiment with the quartz plate and during the 
experiments with the three best tourmaline 
plates. The observations made are plotted in 
Fig. 5, where relative power is plotted against the 
plate voltage on the pentode. Curve 1 represents 
the output of the pentode, quartz controlled. At 
the voltages shown, the tourmaline plates stopped 
oscillating. They would start again on lowering 
the voltage but always stopped on increasing it 
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Fic. 5. Comparison of power output of quartz and 
tourmaline. Curve I represents the quartz. The other 
curves represent the three best tourmaline plates. 


near the values shown. It is realized that a strict 
comparison would require the use of an X-cut 
quartz plate, i.e., one in which the applied field is 
perpendicular to the electric axis as is the case 
with the tourmaline plates. This was not avail- 
able when the experiment was performed. 

In general, one fails to see from these experi- 
ments, wherein tourmaline has any advantages 
over quartz. The cost and scarcity of the material 
are serious objections and the smaller size of the 
crystals make it more difficult to handle. For the 
stabilization of ultra-high frequencies, it still ap- 
pears better to use quartz plates of lower fre- 
quency and triple or quadruple to obtain the 
short wave-length radiations desired. 
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The Theory of Refraction Shooting 


Morris Muskat, Gulf Research Laboratory, Pittsburgh, Pa. 
(Received October 3, 1932) 


Considerations of geometrical optics indicate that the 
limiting rays travelling along the interface between two 
homogeneous elastic media should carry inappreciable 
amounts of energy. Nevertheless, the “first arrival” waves 
in refraction shooting processes, which give linear time 
distance curves, are usually interpreted as being due to 
such limiting rays travelling along the interface with the 
velocity of the lower medium. To clear up the situation 
an analysis of the problem has been made both from the 
points of view of geometrical optics and wave theory. 
It is pointed out that although the assumed “refraction” 
paths are minimal time paths in the sense of Fermat's 
principle the applicability of the latter along a surface of 
discontinuity may be open to question. On the other 
hand, it is proved that, geometrically, they are necessarily 
real, since they are the only type of path that can give the 
observed linear time distance curves. The wave theory 
treatment of the problem for fluid media modified slightly 
from the operator method analysis previously given by 


Jeffreys, is then presented to show that the wave theory 
actually gives waves whose geometrical interpretations are 
exactly those of the refraction shooting process. Their 
amplitudes vary inversely as the square of the distance 
from the source and they give vertical displacements of 


the same order as those of the directly reflected waves.. 


The analysis is extended to the case of general elastic 
media and it is shown that four types of “refracted” 
waves will be produced upon the incidence at the interface 
of either a longitudinal or transverse wave pulse. Two of 
the waves will be recorded as longitudinal waves and the 
other two as transverse. One of each pair effectively 
travels in the second medium with the longitudinal velocity 
of the medium and the other two travel with the transverse 
velocity of the lower medium. The velocities and acceler- 
ations produced by these “‘refracted’’ waves will vary over 
the pulse thickness exactly as the displacements and 
velocities in the incident pulse. 


T is an established practice in applied geo- 
physics to assume that the refraction paths of 
the seismic waves in a “refraction shooting” 
process have a form equivalent to that indicated 
in Fig. 1. 


RECEIVERS 


Fic. 1. Diagrammatic of “refraction” paths 
usually assumed in geophysical analyses. 


meowm (2) 


The feature of these paths which gives them 
the role of refraction paths, and the main concern 
of the present discussion, is, that according to 
such a representation the cone of rays incident 
from a source S on an interface, CD, enters the 
second medium, travels with the velocity of the 
second medium along the interface, and then 
emerges at various points of the interface at the 
angle of total reflection to be finally recorded by 
the seismographs at the surface SR. The reason 


why these paths—apparently satisfying the laws 
of geometrical optics—need discussion, is that 
from an energy point of view they would be ex- 
pected to carry inappreciable, if any, amounts 
of energy. In fact, if medium (2) is taken as 
strictly homogeneous no ray entering it from 
medium (1), even including the limiting one at 
total reflection, can ever again re-enter medium 
(1). Thus a cone of rays impinging at any point 
of CD is, in general, partially reflected and par- 
tially refracted and except for incidence upon 
other interfaces these components suffer no fur- 
ther changes in their paths. In the following, it is 
proposed to review more critically the whole 
situation and to show that the detailed considera- 
tion of the wave characteristics of the problem 
explains the apparent discrepancies. 


THE MINIMAL TIME CHARACTER OF THE 
REFRACTION PATHS 


It is well, however, first to look more closely 
at the practical problem which apparently re- 
quires the assumption of such paths as are drawn 
in Fig. 1. In this regard it is to be remembered 
that these paths are specifically intended to carry 


th 
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the ‘first arrival’ impulses from the sources S 
to the receivers at R. That is, they are the paths 
along which the time of transmission of the 
impulse is a minimum. That they satisfy this 
requirement is readily shown by the following 
elementary calculation. Denoting SR by x, and 
the depth of CD below SR by h, the time of 
travel along a path incident on CD at 4@, is 
evidently: 


t=x/v2+2h{sec 0/v,—tan 0/v2}, (1) 


where, as usual, it is assumed that 22>. This 
has a minimum, 


tmin. = X/V2+2h cos 0,/v1 (1a) 
at the angle of total reflection 
(1b) 


Further, tmin., as given in Eq. (1a), is smaller 
than ¢, which corresponds to the directly reflected 
ray and is given by: 


(x? +4h?)!/?/9, (1c) 
since: 
min. = (1/01?)(x cos 0,—2h sin 6,)?. (2) 


In the sense that Eq. (1a) shows that the 
refraction paths of Fig. 1 represent minimal 
time paths, they satisfy the requirements of 
Fermat’s principle. However, one may question 
the significance of this since the paths in which 
the sections UW are traversed in the first me- 
dium are certainly ‘‘neighboring paths’ in the 
sense of Fermat’s principle, and yet require times 
greater than tnin. by finite amounts. In fact, the 
applicability of Fermat’s principle for paths par- 
allel to surfaces of discontinuity does not seem 
to have been established as yet by any detailed 
analysis and a minimization on one side of such 
a surface may be open to question. But even 
granting the validity of the process used in deriv- 
ing Eq. (1a), its correctness does not of itself 
imply that appreciable amounts of energy will 
actually travel along these geometrical paths. 


THE IMPLICATIONS OF A LINEAR TIME DISTANCE 
CURVE 


In view of this lack of conclusiveness of these 
geometrical path arguments, with regard to the 
reality of the supposed refraction shooting paths 
it is of importance to see how necessary they are 


from an empirical point of view. And here the 
implications are more definite, although the 
totality of the significant empirical data states 
the single fact that in many observed cases the 
time of first arrival at the receivers R varies 
linearly with the distance between the receivers 
and the source, i.e., 


t=ax+b. (3) 


Of course, in many cases there is a curvature to 
the ‘‘time distance” curves but of importance 
here is the fact that there are definitely some 
whose linearity is exact within the experimental 
errors. It should be mentioned, too, that relations 
as Eq. (3) do not hold for all values of x. For small 
values of x the direct wave from S to R will come 
in first, and only after x 
will a wave governed by Eq. (1a) come in sooner. 

The question then is: what, if any, are the 
implications of Eq. (3) with regard to the paths 
along which the “‘first arrivals’’ travel? To answer 
this question it is convenient to set up a general 
relation for the time required for a wave to 
traverse any general type of path. 


>. 


wa 


Fic. 2. Path element of a general seismic wave path. 


Thus let us suppose that the wave reaching 
the distance x, from its shot point, at the time ¢ 
has left the shot point at the angle i) with the 
vertical; and let its angle with the vertical along 
any element of its path be 7. If the velocity 
(variable only with the depth) at this element is 
v, which we may leave perfectly arbitrary, we 
must have by the generalized law of refraction: 


sin 7/v=const.= sin i9/%)=1/a. (4) 


Considering the path as given by the curve 
y=y(x), it is evident that sin i=(1+~’)-"? so 
that together with Eq. (4) 


(4a) 
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The total time and distance for the path are 
therefore given by: 
t= f(ds/v)=af (dy/v(a? —v*)"?) (5) 
x= (vdy/(a? —v*)"*) (6) 
so that 
t=x/a+(1/a) —v*)"?/0 


Since by assumption v depends only on y and 
moreover increases with y, the path is evidently 
symmetrical about the vertical through the point 
of greatest depth along the path; this is given 
by the value of y at which 


v(y)=a. (7) 


(Sa) 


The integral in (5a) can therefore be written 
as a definite integral from 0 to ymax. as: 


x 2 
t= ———— dy. (5b) 
a avo v 
Hence: 
d Ymax dy 
— (ta)= 14200’ 
dx 
and finally 
dt/dx=1/a. (8) 


Hence if ¢ varies linearly with x as in Eq. (3), 
a must be the same for the paths terminating 
at the different values of x. This, however, is 
impossible if the path is to have a single point 
minimum, as implied by Eq. (5b), and as one 
would conclude from the observation that once 
a and hence ip is chosen the complete path with 
its lateral extent, x, is immediately determined. 
Eq. (8) nevertheless demands that if the time 
distance curve is linear, as in Eq. (3), a given 
angle of entry ig must give rise to a number of paths. 
Hence, since we must still maintain the require- 
ment of symmetry of the paths about the mid- 
way plane, it must be concluded that these 
various paths arising from the single angle 7» or 
single a are extended horizontally at the depth 
of their maximum penetration, just as the refrac- 
tion paths of Fig. 1 are supposed to be. Or, 
stated more explicitly: No set of seismic paths 
suffering refraction due to a velocity distribution 
increasing (Continuously or discontinuously) with 
depth, can give rise to linear time distance curves 


unless the paths have a common maximum pene- 
tration and continue horizontally, for variable 
distances at the maximum penetration. 

Since as already mentioned, linear time dis- 
tance curves are actually observed empirically, 
the above discussion may be considered as a 
proof by means of geometrical optics of the 
reality of the refraction paths of Fig. 1. These 
paths, then, do not represent simply a fortuitous 
assumption, but rather the necessary implication, 
from the point of view of geometrical optics, of 
the observed linear time distance curves. ; 

It appears, then, that although it is difficult 
to see from the elementary consideration of geo- 


‘metrical optics principles, how the refraction 


paths can be of significance with respect to the 
transfer of energy, still the empirically estab- 
lished linear time distance curves do demand the 
postulation of these paths. To resolve this diffi- 
culty it is evident that one must make a close 
analysis of the details of the whole process from 
the wave point of view. 

Before entering upon such a discussion, how- 
ever, it may be of interest to suggest a physical 
interpretation of the refraction paths of Fig. 1 
from a somewhat different point of view than 
has been heretofore considered. 


A 


Fic. 3. Plane wave representation of the ‘‘refraction shoot- 
ing’’ process. 


Thus, as indicated in the diagram, a wave 
travelling with the velocity v at an angle 6 with 
the normal has a horizontal component velocity 
of vsin 6, but at the same time an apparent 


¥ 


~. 
~ 
\ As. 
\ k 
/ 
\ 
\ / 
\ / 
\ / 
\ 
\ J 
\ / i 
/ 
\ 


CY 


REFRACTION SHOOTING 17 


velocity as measured, for example, by the time 
required for an impulse, on the same wave front, 
to reach B after it has reached A, which is given 
by: 

v,=v,/sin 0. (9) 


If then we consider a single plane wave pulse 
passing a point as O and travelling towards the 
interface to be reflected and observed at A and 
B, the effective velocity of propagation between 
A and B will be given by Eq. (9). Now it may 
be shown that the time required for a plane 
wave impulse passing O to reach a point as B 
by a direct reflection at the interface is a maxi- 
mum for the regular angle of reflection, i.e., for 
the ray which would meet the interface midway 
between O and B, and is evidently a minimum 
for a vanishing angle. In any case, if only the 
times for pulses travelling at the same angle are 
compared, the time distance curves found at A 
and B will always be linear. And if bundles of 
pulses moving at different angles are available, 
those with the smaller angles should come in first. 
However, since the reflected amplitudes will be 
greater when the angles equal or exceed that of 
total reflection the first arrival pulse of appre- 
ciable amplitude may well be that incident at 
the interface at the angle of total reflection. 

Since for the angle of total reflection: 


0,=v2 (9a) 


the first arrivals of appreciable amplitude would 
have the apparent velocity of the second medium. 
Further, as already mentioned, they would give 
a linear time distance curve. _ 

In this way the main features of the refraction 
paths may be developed without any recourse 
to the waves actually travelling in the second 
medium. The refraction paths are, then really 
all reflection paths except that observations at 
successive stations are to be interpreted as relat- 
ing to the same wave front. In this sense this 
point of view represents another extreme idealiza- 
tion. Thus whereas previously all the paths have 
been represented by concentrated line trajec- 
tories passing from the sending to the observing 
points, the above picture is based on the concept 
of infinitely extended plane wave fronts, different 
parts of which may be recorded at various sta- 
tions at correspondingly different times. 


Although this point of view does avoid the 
question of energy transfer in the ray travelling 
along the lower surface of the interface, it, too, 
is evidently greatly idealized and cannot give 
a complete description of the phenomena in- 
volved. In particular, the infinite plane waves 
are simply impossible in a finite system and 
certainly near the wave source the assumption is 
quite unjustified. On the other hand, the above 
point of view does serve to bring out the fact 
that all that we really observe at the surface are 
the superpositions of a number of individual 
wavelets in the form of a wave front, rather than 
a series of individually projected rays which can 
be separately traced back to the original 
source. And in fact, the wave theory analysis 
which we shall now present does give just this 
type of a representation. 


WAVE THEORY ANALYsIs FOR LONGITUDINAL 
WAVES 

The wave theory analysis of the problem of 
refraction shooting has really already been given 
by Jeffreys' for the separate cases where only 
longitudinal or transverse waves are considered, 
although Jeffreys’ interest in this question arose 
from problems of earth quake propagation. In 
the following, to show the method of analysis 
that will be used later for the more general case, 
Jeffreys’ discussion of the problem for compres- 
sional waves will be essentially repeated although 
his operator method will be replaced by the 
somewhat more familiar differential equation 
analysis. Further, since we are primarily inter- 
ested in the effect of the interface we shall from 
the outset consider the simplified problem of a 
single interface separating two infinite media. 
Jeffreys treats the case where the upper medium 
is bounded but then neglects the waves reflected 
from the upper free surface. On the other hand, 
we shall analyze the reflected wave in somewhat 
more detail than is given in Jeffreys’ work, and 
finally, in the last section we shall treat the case 
of general elastic media which reveals the exis- 
tence of several types of ‘‘refracted’’ waves. 

In the present section, however, our problem 
will be to find the waves resulting from a spher- 
ical pulse emitted at a source in a homogeneous, 


'H. Jeffreys, Camb. Phil. Soc. Proc. 23, 472 (1926); 
Gerl. Beitr. z. Geophys. 30, 336 (1931). 
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non-dispersive, and non-aborbing fluid medium 
at a height above the interface, z= 0. Because 
of the symmetry the cylindrical coordinates (9, z) 
will be used. The densities of the media will be 
denoted by y:, v2, and the velocities by (;, v2) 
where 

Since only compressional waves will be treated 
for the present, the problem may be described 
by the potential where: 


(10) 


and the displacement components along z and p 
are given by: 


w=db/dz; q=db/dp. (10a)? 


Eq. (10) implies that the problem is in a steady 
state and composed of waves of the single fre- 
quency w. The problem of the pulse will be 


derived by a superposition of the effects of the | 


individual frequencies. 
Fundamental solutions of Eq. (10) are of the 


form: 
P= et** Jo(Ep), 


a?— 2 


where (11) 
and these must be combined so as to give a 
source at (0, 2) and also to satisfy the boundary 
conditions at z = 0. To obtain the effect of the 
source it is more convenient, however, to simply 
add on directly the function 


e~ + (h —z)?}/? 


which satisfies Eq. (10) and represents a unit 
source at (0, 2), emitting compressional waves of 
frequency w and velocity v. But since for the 
remainder of the solution, functions of the form 
of Eq. (11) will be used. Eq. (12) will also be 
expressed in this way. Thus from well-known 
results of the theory of Bessel functions we may 
write: for z<h 


Jo( Ep) EdE 
Po= if ’ 
0 


a 


(12)8 


(12a) 


2 has the additional physical significance that for a 
harmonic wave it is proportional to the dilatation of the 
medium. For a general wave form the dilatation is pro- 
portional to 0*@/dt*. Further, its time derivative is a 
velocity potential and is the function Jeffreys uses. 

3 For w <0, all the succeeding formulae are to be changed 
to their complex conjugate values. 


where a@ is again given by Eq. (11), with v=2). 
The complete solution to the problem may 
now be expressed as: 


z>0 

0 
2<0 


(13) 


4 


where a?= #—*/v2’, and A and B are arbitrary 
functions to be chosen in such a way that the 
normal displacements and pressures are con- 
tinuous at the interface z=0, i.e., such that: 


W1= 0D, /dZ= A’2/dz 
atz=0. (13a) 
i= 


If Eqs. (13a) should hold, A and B must be 
given by: 

B= 2te-*/(ak+a); 

A = (fe~*"/a) (ak —a@)/(ak+@); 


Hence we have finally: 


k= 2/71. 


ak—a@ ) 
70 k+a 


aa 


2 Jo(Ep)dé. 
0 ak+a 


To interpret Eqs. (14) we shall again follow the 


(14)* 


*It may be of interest to observe that although equations 
exactly equivalent to Eq. (14) occur and have been dis- 
cussed in detail in the problem of the transmission of radio 
waves over the earth’s surface (cf. Riemann-Weber, Die 
Differential und Integralgleichungen der Mechanik und 
Physik \1, p. 541) the fact that these analyses have not 
led to “refraction paths’’ such as we shall find below, is in 
no way contradictory to the analysis to be given here. For 
in the problem of radio-wave transmission the complex 
velocity in the earth is in absolute value Jess than the 
velocity in the incident medium, air. Hence, one would not 
expect, even geometrically, total reflection phenomena and 
the associated refraction paths which arise in the seismic 
problem. And in fact the approximations made in evaluat- 
ing the integrals of Eq. (14) for the radio-wave problem are 
explicitly based upon the assumption that the effective 
velocity contrast is just the reverse of what it is in the 
seismic case. In this way, together with the assumption of a 
complex effective earth velocity, our refraction paths are 
not only automatically eliminated, but at the same time the 
Sommerfeld surface waves are introduced as the predom- 
inant feature of the phenomena. 


) * 
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method of Jeffreys and of Sommerfeld.‘ Thus, 
first, it is convenient to separate out the second 
term of #; since it alone represents the reflected 
wave or the effect of the second medium. At the 
same time the change of variable: &=(w/v,)n; 

= (w/v1)(n?—1)"?; @= (w/01) (9? may 
be introduced. We have then: 


n ak—a@ wp 
en atht2) J, dn. (15) 
aak+a V1 
Evidently the integrand has branch points at 
a=0, a=0. To define the behavior of the inte- 
grand as these points are crossed we shall take a 
and @ as positive when 7>1 and make a cut in 
the 7 plane so that » moves around n=1 and 
n=1,/V2 on the positive imaginary side. Eq. (15) 
may now be broken up into 2 integrals by the 
substitution of the relation: 


2Jo= TT +11), 


where HH) and H)® are Hankel integrals, the 
former vanishing on an infinite semicircle in the 
upper half of the complex plane and the latter 
when the semicircle is in the lower half of the 
complex plane. Since the integrand is regular in 
the upper right quadrant of the complex plane 
and its branch points lie below the cut along the 
real axis, Eq. (15) may be replaced by: 


w n ak—@ 
——— e-atht2) FT (1) (= =) dn 
aak+a 
n ak—@ 
en atht2) 7,2) dn 
20)? J-tec a ak +a 
n ak—@ 


wpn 
e~atht 2) JJ (2) (= dn, 


2v;? aak+a V1 


where the circuits in the last integral are infinite 
loops about the branch points, and which close 
at the infinite quarter circle. 
Since Ho (iz)=—H)®(—iz), the first two 
integrals above cancel and we are left with 
n ak—@ 
——_ 
aak+a 


wpn 
enathtz) (16) 
V1 


Considering now the loop around the branch 


point at a=0, the contribution of the small con- 


* A. Sommerfeld, Ann. d. Physik 28, 665 (1909). 
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tour about the branch point itself vanishes as the 
loop is shrunk and one has simply the effect of 
a changing sign as the branch point is encircled. 
Thus we may set 


n=1-—10:0>0. (17) 


Now since /J)®(z) vanishes exponentially for 
large negative imaginary parts of its argument, 
the integrand is small except when @ is small. 
We can therefore set: 


a 


B= (1 —v,"/v27)"/2, (17a) 


As the branch point, a=0, is circled and the 
contour returns to the infinite quadrant the 
sign of a is changed. Expanding the term 
(ak—a)/(ak+a) in powers of a/@ we have to 
second order terms: 


n 


vy 0a 
2a 
{cosh a(h+z)+— k sinh 
a 


Replacing now, for p large, 7) by the first 
term of its asymptotic expansion, and setting: 
m=a(h+z)/6'?, we may rewrite ®, as 


xX {cosh my+[2km/a(h+z) sinh my}dy 


1 mkv, -(18) 

p awp(h-+2z) 

p 


where the term in 1/p? has been dropped since 
an equivalent approximation has already been 
made in the substitution for //)®. 


5It is to be noted that this approximation and the 
corresponding one for @ in Eq. (19) really imply that 8? 
is much larger than the important values of @ and this 
means, physically, that there is a large velocity contrast 
between the two media. The general features of the 
problem, however, will not be changed by this implied 
assumption. It may also be mentioned that Eqs. (17) 
(17a) imply that the loops are to be drawn parallel to the 
imaginary axis. 
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In a similar manner the contribution #,, due 
to the branch point at @=0 may be found. Thus 
setting: 


V1 tw 
n=——10; B; 
v v 
2 1 (19) 
V1 1/2 
es 
V1 Ve 4 


a has the same value on the two sides of the loop 


1 


—ie a/v, 


while @ changes sign. The final result is: 


©,,= — (20) 


THE REFLECTION OF PULSES 


The physical meaning of the two terms #, 
and #,, may be seen more clearly by forming 
pulses of these harmonic waves to correspond to 
an initial pulse formed from #9. Thus we may 
define the incident pulse by the Fourier integral: 


(21) 
+ (2—h)?)"? ic 0 
which represents a pulse defined by 
4+ =- (21a) 
(p? + (h—sz)?)'? V1 V1 Vv) 


=( : otherwise. 


The addition of —ic to the limits of the integral 
in Eq. (21) has been made to take care of con- 
vergence requirements in the integrals for ®, 


that follow. 
Performing now the same operation on #,, and 


®, as carried out on ®o, we obtain 


+cao—ie 
dw J el y, (22) 
2rp 
—iv;? +co—ic dw a/v, 
_w-ic w Yo 


By comparing ®,, with it follows that: 
—(1/p)fl — t+ | 

Recalling the fact that in deriving Eq. (18) 
we assumed that p is large it is evident that 
p+(h+z)?/2p is really an approximation to 
(p?+(h+:)?)"*. To this approximation therefore, 
#, may be finally rewritten as: 


,,= —(1/p)fL—t+(e? + 
205 


=0: otherwise. (24)° 


By comparing the argument of f with Eq. (1c) 
it becomes clear that ®,, represents the wave 
directly reflected from the interface according to 
the laws of geometrical optics. Thus ®,, is zero 
until the time that would be required for a ray 


° @,, as given by Eq. (24) corresponds to the amplitude 
of a wave fotally reflected from an interface. This is 
.evidently due to the implicit assumption, already men- 
tioned above, that the velocity contrast across the inter- 
face is high. 


to pass from the front of the pulse (at t=0) 
along the ordinary reflection path to the point 
(p,z). Furthermore, this reflected wave is a 
pulse of exactly the same form as the incident 
pulse. Hence , behaves as if it were emitted 
by the negative image of the real source in the 
interface, so that to the approximation that 
p~(p?+(h+2)*)'? the radial displacements and 
velocities of &) and ®,, are equal and opposite. 

The vertical displacements and velocities of 
this reflected wave are, however, greater than for 
the direct wave ®o. This will be clear from the fol- 
lowing list of the various displacement compo- 
nents in which coefficients as 1/(p?+(h-+:)?)'/ 
have been replaced by 1/p and only the terms 
in lowest power of 1/p have been retained: 


(h+2) } 
= —t+ ). 
vip" V1 
(25) 
1 
f —t+ ). 
vip V1 


| 


‘1 
8) 
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whereas: 


wvo= 


vip? 


1 (p?+(h —z)?)12 
go= ). 
J 


(p?+(h } 


V1 
(26) 


It may be finally noted from these that the 
radial components decrease only as 1/p whereas 
the vertical components decrease as 1/p*?. And 
this is to be expected since for large values of p 
the reflected wave must approach in form a 
longitudinal cylindrical pulse in which the radial 
components naturally predominate. 

Returning now to ®,,, we see clearly that since 
it has the term 1/w as a factor in the amplitudes 
for the individual harmonic waves, it must repre- 
sent a wave distorted with respect to the initial 
pulse. Its pulse character, however, becomes im- 
mediately evident when the line contour of Eq. 
(23) is completed by an infinite semicircle in the 
lower half plane. For then, by Cauchy’s theorem, 
it follows at once that: 


or if 


if (27) 


(27a) 


Recalling now Eqs. (1a) and (17a), and noting 
the fact that a/v, is simply the time required for 
the wave pulse to travel a distance equal to its 
thickness it is clear that %,, vanishes until a 
time required for the pulse to traverse the refrac- 
tion path of Fig. 1. In other words, ®,, represents 
a wave that has apparently travelled along the refrac- 
tion path of Fig. 1. The reality of these paths may 
therefore beconsidered as conclusively established. 

Although this has been the essential aim of the 
above analysis it is of interest to note some of the 
features of this refracted wave. Thus as already 
suggested #,. itself is distorted with respect to 
®). However the velocity potential and displace- 
ments the same shape as Thus: 


OP, 
ot kp 23? p? Vo V1 


— 20 B(h+3) 
(28 


2v;? 
qr,= 
Bp? Pp V1 


all with the restriction 


Comparing these with the corresponding val- 
ues for po in Eqs. (21a) and (26) it will be seen 
that if @) has a sharp front so will the displace- 
ments of the refracted waves. On the other hand, 
if the displacements in the incident pulse rise 
sharply as the pulse is entered, the velocity com- 
ponents of the “refracted” pulse will rise sharply. 

As to order of magnitude, Eqs. (28) show that 
both the radial and vertical components of dis- 
placement and velocity fall off as 1/p*. Hence, 
whereas in the incident or directly reflected pulse 
the radial components bear an increasing ratio 
to the vertical components as p increases, in the 
refracted wave they bear a constant ratio. Thus: 


p/(h+z), 
qr,/Wr,= 01/V28= tan 6, 


while (29) 


where @, is the angle of total reflection. g,, and 
w,, are therefore actually of the same order of 
magnitude. Incidentally, this last relation shows 
explicitly that the refracted waves travel in the 
first medium at the angle of total reflection just as 
is shown in Fig. 1. 

The relative magnitudes of the refracted and 
directly reflected waves are indicated by the 
ratios: 


qr, 28% | 
Wry (31) 


RveB(h+z) 


It follows that the reflected wave will give 
the major radial effect at large distances, but 
that both types of waves contribute appreciably 
to the vertical components. In fact, taking the 
function f and f’ to be of the same order of 
magnitude, k= y2/yi~1, v1/v2~1/2, and s=h it 
follows that w,,/w,,~0.580,/h. Hence, in this 
case, unless the source of the incident wave is 
distant from the interface by more than half of 
the wave velocity in the upper medium, the 
vertical displacements due to the refracted wave 
will even exceed those of the directly reflected 
wave. Of course, it will be understood that these 
comparisons between the various waves do not 
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refer to the same instant since at the very begin- 
ning of this discussion it was pointed out that 
for large distances the refracted wave would 
arrive at a given point sooner than the reflected 
wave. Rather, they refer to the effects of the 
waves at the same points at whatever times the 
individual waves may arrive there. 

It may be concluded, therefore, that the wave 
theory not only confirms in every detail the 
reality of the paths of Fig. 1 but it shows further 
that the refracted waves travelling along these 
paths have vertical displacements of the same 
order of magnitude as those of the directly 
reflected waves, and hence will have a vertical 
energy flux of the same order as that of the 
directly reflected waves. 

The following remarks, however, may be 
added. First, the analysis presented here has 
treated the case where the initial wave source 
is embedded in an infinite medium whereas the 
actual problem involves a free surface—the 
earth’s surface. It is true that the detailed results 
would be modified but the essential features and 
existence of the refracted wave would evidently 
remain unaltered. In fact, the free surface can be 
readily taken care of, as Jeffreys has done, by 
adding to the real wave source its negative image 
in the free surface. The main new feature intro- 
duced in this way is the set of pulses reflected 
from this free surface but which have no bearing 
on the problem of ‘‘refraction shooting.” 

Another point that may be mentioned is that 
the approximation made in the above analysis 
that p has to be large really refers to the require- 
ment that wp/v; or 27p/X, is of the order of 1 
or greater, where \ is the wave-length in the 
upper medium. Since p is always large when the 
whole phenomenon refers to earthquake dis- 
turbances the approximation will always be valid 
there. But in geophysical work where p/v; may 
be as small as 1/2, the approximation may seem 
questionable. However, in practical cases the 
predominant frequencies are usually of the order 
of 50 cycles per second or greater and then the 
product wp/v, will again be large and the approxi- 
mation will again be justified. 


GENERAL ELAsTIC MEDIA 


Somewhat more serious, however, than these 
limitations of the above analysis is the fact that 
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it has referred to fluid media, in which the trans-: 


verse waves are automatically eliminated. Again 
the reason for this simplification has been the 
fact that the integrals in the general case that 
correspond to Eq. (14) are more difficult to 
evaluate completely. But since the analysis that 
can be carried through brings out some features 
that might not be anticipated, it will be briefly 
outlined, although in one respect the solution 
will still be incomplete. 

For the general case where both longitudinal 
and transverse waves are excited at the interface 
one must add to the previous function , another 
function x, to take care of the transverse wave 
components. 

The radial and vertical displacements for a 
harmonic wave system of frequency, w, are then 
given by: 


r (32) 


where v and 3 are the longitudinal and transverse 
wave velocities respectively and a common factor 
e*' has been dropped throughout. Denoting by 
subscripts 1, 2, the media 1, 2 (separated by the 
plane z=0, cf. Fig. 4) the wave system resulting 
from a longitudinal wave source at (0, #) in 
medium 1, may be represented as: 


+ 


(33) 


A 
0 


Bee *Jo(Ep)dé, 
0 


B= (33a) 


| 
Op Ob Ax | 
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It will be supposed that 7,;<v;<i2<v2. continuity of the displacements and stresses at 
Expressing the first term in #; as in Eq. (12a) | the interface, the following equations are ob- 
and applying the boundary conditions of the | tained: 


—A; +8:B, +A2 +B2B2 
—a,A,;+#2B, a fem", 


(34) 
where 6= €= 22? wo, being the rigidity moduli of the media. 
The determinantal denominator of the solutions for A,, B,, As, Be is: 

A= ki? — €1) }? — +22(6—1) +8201) \ (34a) 

AA yaye%"/E= ky? — (b€2— €1) +282 +28 (6 —1) — (34b) 

X [beg —2#(6—1) —ki?} — 1)? — 
AB 2a; {| 2a2B2(5 — 1) — €1) — (Gee — €1) (See — 28) (34c) 
with similar expressions for Ao, Be. | the integrands of the integrals in Eq. (33) will 


As the above equations make it possible to have branch points in the complex plane at: 
write out explicit expressions for P2, x1, x2 
the problem is now formally solved. ,A useful ae 
although the complexity of A makes it difficult 
A to be certain as to whether or not there are any 
poles for positive values of the real part of & In 
this respect, the present analysis will be incom- 
plete. Nevertheless the contributions of the 
branch points and the general nature of the 
\) contributions of the poles, if there are any, can 
be readily derived. 
¥Y Thus it easily follows as before that the in- 
tegral in ®; can be reduced to 


> 
$,= —(w/20:) £Ai(n) 
(2) (wpn/vi)dn tT, (35) 
“Y where the substitution £=w7/v; has been made 


in Eq. (35), so that a;=w/v,(y?—1)"*. The cir- 
cuit for the integral is composed of loops which 
come from the lower right infinite quarter circle, 
encircle the branch points and return to the 
quarter circle. 7, represents the contributions 
due to the poles. 

Now when the loop encloses a particular 
branch point no=21/v0, Ai(n) may be written as: 


Fic. 4. Coordinate system for wave theory analysis. 


discussion, however, requires the evaluation of 
the integrals. To carry out this evaluation one A n(E+ Fao)e~™" 
may still use the method of contour integration i(n) = 2 1/2( BY BY 
(4? —1)"?(E’ + F’ao) 
used in the simpler case of pure compressional 
waves. Thus it is quite clear from Eq. (34a) that | where ao= (n?—7,*)'”, and is to be considered as 


(35a) 


| 

1 | 
) 

| 
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small in comparison with E, F, E’, F’ which 
may be taken as constants over each loop. For 
the loop about , 7 may be replaced by mo—i0 
and then a» takes the form ap=(2v;/v9)'/*e-*"/4, 
Further, as the branch point np is encircled E, F, 
E’, F’, are approximately constant while ap 


®,, (E/ 22/200 


where N is the numerator of- A,(n), except for 
e-™'", and mo represents the position of a pole A. 
The upper or lower signs are to be used as 7 
is above or below the real axis. The F, E, F’, EF’ 
for the different ®, are of course different and 
may be found from Eqs. (34a) and (34b). 

On recalling the discussion of Eqs. (18), the 
following geometrical interpretation of ,---7 
may be readily verified. The wave ,, will arrive 
at (p, z), according to Eq. (36a), at a time: 


U,,)= (37) 


after the incident wave has been started at (0, /). 
It is, however, a wave damped with respect to 
its distance from the interface and the distance 
of the source from the interface. It must therefore 
be excited at and travel along the interface. It 
is apparently excited initially by a direct longi- 
tudinal wave to the interface, travels along the 
upper part of the interface as a transverse wave 
and excites vertically projected and damped 
longitudinal waves. This history of #, would 
imply that its time of arrival would be 


rather than that given by Eq. (37). But since it 
has already been assumed that p is large and the 
velocity contrasts are high it is seen that Eq. 
(37) may well be taken as the approximation to 
Eq. (37a) in which the term [(4+2)/p ]d;/v; has 
been dropped because it is small compared to 
unity. 

®, is simply the directly reflected longitudinal 
wave and gives nothing new. Its time of arrival is: 


t= (®,,) = (p? + (38) 


(37a) 


changes sign.’ Putting in specifically 
going to large p, and replacing JJ, by its 
asymptotic expansion, we find that to the lowest 
power of 1/p: 


?,= 71, (36) 


where: 

(36a)8 
(36b) 
(36c) 
(36e) 
(36f) 


#,, however, is also a new wave. Its time 
of arrival corresponds to paths in the upper 
medium, at the angle of total reflection for the 
refracted transverse waves, which are joined by 
a path on the under side of the interface trav- 
ersed with the transverse velocity of the lower 
medium. ,., itself, of course refers to the wave 
that is observed at (p, z). Its origin is apparently 
the cone in the incident wave that travels at the 
angle of total reflection for the refracted trans- 
verse waves. The next stage in its history (as 
suggested by its time of arrival and angle of 
emergence into the second medium) corresponds 
to the “totally refracted’”’ transverse wave ex- 
cited by the incident wave and which travels 
along the under side of the interface. And this 
finally is partially transformed into the longi- 
tudinal waves #,. reemerging from the interface 
at the same angle as the original exciting wave. 
That this angle of emergence really corresponds 
to that of total reflection, 6;,, between the longi- 
tudinal waves and refracted transverse waves, 
readily follows from the fact that by Eq. (36c) 


r,/Wr,= 01/021 = tan Ou. 


7It may also be noted that EZ, F, E’, F’ are all pro- 
portional to w* so that except for e~*"”, A,(n) is independent 
of w. 

8 When w has a negative real part, the conjugates of the 
above , should be used, except that factors as e~*” must 
be replaced by e“”. These changes, though fairly obvious 
from general physical considerations may be derived 
rigorously by beginning with an incident wave as Eq. (12) 
with w negative, putting the indentation in Eq. (12a) 
above the real axis and evaluating the loops which are 
drawn in the first quadrant. 


st 


REFRACTION 


The time of arrival of this wave is: 
t(®,,) = p/de+L(h+2) 


And again, ®,, is the ‘refracted’ wave that 
appeared in the simpler case of the purely com- 
pressional wave system. It behaves as if it travels 
in the second medium just below the interface 
with the longitudinal velocity in the second 
medium and travels in the first medium at the 
angle of total reflection between incident and 
refracted longitudinal waves, 6). Its time of 
arrival is given by: 


t(?,,)= p/ve+[(h+s)01 


(39) 


(40) 


71, if it exists, is a real surface wave provided 
the imaginary part of y is small. For then it will 
decrease essentially as 1/p'/? which is the mode 
of decrease for a surface wave. It travels along 
the upper side of the interface with an effective 
radial velocity equal to the longitudinal velocity 
in the first medium divided by the real part of no. 
A diagrammatic representation of these waves 
is given in Fig. 5.° 

These waves all refer to the dilatational or 


Xr, = (iw E / E” 1) ( E's+0,F’ /w) /21p—(wh/ ry) 


T= N(no)/A’ (no) (oh! (ne? —1)U2—(w 2/04) (ne? i(wpn 


where in the expressions for x, the E and F are 
to be taken from the coefficient of a; in Eq. (34c), 
so that E, F, ~w'/zv;' and E’, F’~ 

The geometrical interpretation of these is the 
following. x,, is a transverse wave that has been 
apparently excited at the point on the interface 
that is just below the source. The time of arrival 
for x,,, according to Eq. (41), simply includes 
the time of transit after excitation, 

U(x) = (e? (42) 

* When synthesized into pulses #,, will have the same 
form as % #,, and %,, pulses, however, will be distorted 
in such a way that the velocities and accelerations of ®,, 


and ®,, will respectively have the same forms as those of 
the displacements and velocities of p. 
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Fic. 5. Geometrical representation of the longitudinal 
waves excited by an incident spherical longitudinal pulse. 
®,,: damped wave; directly reflected wave; ®,,, ®,,: 
“refracted” waves. 01, 0 are angles of total reflection. 
v, v are longitudinal and transverse wave velocities. 


longitudinal components of the disturbances ex- 
cited by the incident longitudinal wave. In addi- 
tion there are of course similar transverse com- 
ponents represented by the functions x. Without 
repeating again the details x; may be represented 
as the sum of the following functions: 


(41) 
(41a) 
(41b) 
(41c) 
(41d) 


but lacks the term //v, that would be required 
for the incident wave to reach the point of ap- 
parent excitation. This, however, seems to be 
due to an implicit approximation in which this 
term, //v,, is neglected, since if it were large the 
exponential damping factor 
would anyway practically extinguish the wave. 
X», appears to be a transverse wave that has 
been excited by a direct longitudinal wave at the 
point (0, p). From Eq. (41a) it would seem that 
it travels upwards from its point of excitation on 
the interface, at an angle with the vertical of 
1/(v,/i2—1)'*. But since, as already 
mentioned it is implicitly assumed in the analysis 
that the velocity contrasts are all large, the 
above angle is the corresponding approximation 
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to a vanishing angle. Its time of arrival is: 


(p?+h?)!/? 2(1 
t(x,,) = + 
V1 v1 
(p?+h?)'/? 
(43) 
V1 V1 


xr, is again a “‘refracted”’ wave. It appears to 
be excited by the incident longitudinal cone’ 
coming at the angle, 6, of total reflection be- 
tween longitudinal and transverse waves; it then 
travels on the lower side of the interface as a 
transverse wave and then reemerges as a trans- 
verse wave at the angle of total reflection be- 
tween transverse waves, 6; Thus its time of 
arrival is given by: 


= (p/V2) +(h/v:)(1 
+(2/0:)(1 


Finally, x,, is another “refracted” wave whose 
exciting cone is incident at the interface at the 
angle of total reflection between longitudinal 
waves, 4;;. Its immediate source is the “totally 
refracted” longitudinal wave which travels on 
the underneath side of the interface and excites 
Xr, aS a transverse wave travelling at the angle 
of total reflection between a transverse and longi- 
tudinal wave, 0. Its time of arrival is: 


= (p/v2) +(h/v1) (1 
+ (2/01)(1 


T is again a surface wave and of no immediate 
interest. 

The waves x, are represented diagrammatically 
in Fig. 6. It should be pointed out, however, that 
only the solid parts of the paths, or those parts 
for which p is large, are definitely indicated by 
the above analysis. For the angle with the ver- 
tical of the direction of propagation of the waves 
xr are given by 0=tan~'w/q (for the waves 
6= tan! g/w) and since the expressions for x, of 
Eqs. (41) to (41d) refer to large p, the corre- 
sponding angles @ will also refer to the regions of 
large p. However, the times of arrival as given 
by Eqs. (42) to (45) indicate what are at least 


(44) 


(45) 


This is the same “wave cone” that excites the 
wave. Similarly, the x,, and ®,, waves are excited, from a 
geometrical point of view, by the same incident ‘wave 
cone,” as may be seen by comparing Figs. 5 and 6. 


MUSKAT 


equivalent geometrical representations of the 
paths which are finally completed by those por- 
tions represented by the analytical, expressions 
for the x,. These dotted portions, therefore, do 
not necessarily represent the analytic continua- 
tions of Eqs. (41) to (41c) for smaller values of p 
but rather the probable preliminary longitudinal 
or transverse waves which finally excited the x, 
waves represented by the solid lines or by the 


V, V, 


Fic. 6. Geometrical representation of the transverse 
waves excited by an incident longitudinal pulse. x,,: 
damped wave: x,,: reflected wave. x,,, x,,: “refracted” 
waves, 911, Ou, 91, are angles of total reflection. v, 
are longitudinal and transverse wave velocities. 


expressions of Eqs. (41) to (41c). The same 
remarks apply, of course, to the paths indicated 
in Fig. 5. 

It may also be of interest to note that the 
paths of Figs. 5 and 6 are actually minimal time 
paths from the point of view of geometrical 
optics. Thus by the method indicated on p. 15, 
it may be readily verified that by minimizing 
general paths of the types of Figs. 5 and 6 not 
only will the angles come out to be the various 
angles of total reflection given by the above 
equations of the wave theory analysis, but the 
minimal times will also coincide with those given 
in the above equations. 

Since the times of arrivals for the above waves 
are different the resultant wave system at the 
point of observation (p, z) will not be a simple 
superposition of all the waves recorded at the 
same instant. Rather, if p is sufficiently large 
and the velocity contrasts are relatively high a 
single longitudinal pulse originating at the point 
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(0, 2) will lead to a series of pulses in the order of 
the times of arrival given in the above equations. 
Since: 


t(®,,) < < < < ta < 
< < < xn) 


the pulses will arrive at a given point in the order 
of the arguments of the above inequality. fa is 
the time for the direct incident pulse and is 
evidently equal to (p?+ (h—z)"/?/v,. 

Of course, after the passage of any particular 
pulse, the disturbances will not suddenly vanish 
so that in an actual observation it may be diffi- 
cult or impossible to distinguish them. But the- 
oretically new impulses should be observed at 
each of the instants indicated by the above values 
of t. The strengths of these impulses will also be 
variable and some may be too weak to be 
observable. These details as well as the shapes 
of the pulses can, if desired, be determined from 
the equations already given."' Here, we are con- 
cerned mainly with the fact that the refraction 
paths are real and that the displacement com- 
ponents are generally of the order of 1/p? and 
hence may under the proper conditions carry 
quite appreciable amounts of energy. 

Finally it may be noted that the above wave 
theory analysis shows that when the incident 
wave is transverse, a similar series of refracted 
rays will arise. Those that really have refraction 
characteristics as ®,., ®,, xX;,. Xr, Of Figs. 5 and 
6, are shown in Fig. 7. The times of first arrival 
are: 


(46) 


ph 
(1 
Vo V1 


h 
U2) =" += (1 (1 L047) 
v 


2 1 1 


= 
t(x2) = p/d2+[(h+2) 
with 


t(1) < t(x1) < (2) < t(x2). (47a) 

" It may be mentioned, however, that when synthesized 
into pulses the real “refraction” waves Xr, and x,,, just as 
the +, waves, will have the same velocity and acceleration 
distributions respectively as the displacement and velocity 
distribution in the incident 4) pulses. 
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It may be observed that if h=< and p is very 
large the above waves have the same times of 
arrival as those for the incidence on the interface 
of a longitudinal wave. Thus we have: 


t(x1) 2 t(®,,); t(x2) 2t(®,,). 


To be certain, therefore, of the nature of the 
incident wave, it would be necessary, if the 
question arises, to determine if the recorded wave 
is itself longitudinal or transverse. Thus if a 
longitudinal wave has been found to arrive at 
the time corresponding to ¢(x,,), Eq. (48) shows 


an 


(48) 
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Fic. 7. Geometrical representation of the “refracted” 
longitudinal (#) and transverse (x) waves excited by an 
incident transverse wave. 011, 91, 9, 9. are angles of total 
reflection. v, 0 are longitudinal and transverse wave veloci- 
ties. 


that it must be the ,; wave resulting from the 
incidence of a transverse wave pulse, rather than 
the x,, wave excited by the incidence of a longi- 
tudinal wave. 

Aside from the significance of these results in 
the theory of refraction shooting it may be of 
interest to observe that they also suggest an 
explanation for the apparent velocity contrast 
that has been reported in the literature’ between 
wave velocities parallel and perpendicular to the 
bedding planes. Thus, rocks or sands which show 
marked ‘“‘bedding planes’”’ may be considered as 
composed of stratified layers each essentially 
isotropic, but varying among themselves in their 
wave velocities. Now the analysis developed 


” B. McCollum and F. A. Snell, Physics 2, 174 (1932). 


s 
O 
- 
| 
| 
f 
0,2) 
Y w/ 
e 
d 
e 4 
e 
i 
g 
it 
s 
e 
n 
e 
e A 
e 
e 
a 
t 


28 MORRIS MUSKAT 


above shows that when a high velocity medium 
bounds one of lower velocity, waves originating 
in the latter will be refracted into the high 
velocity medium and effectively travel in it along 
an interface. This will evidently remain true if 
the high velocity medium is bounded again by a 
low velocity layer unless the wave-length of the 
wave is much larger than the thickness of the 
high velocity stratum. And since in explosion 
waves many frequencies are present it does not 
seem unreasonable that at least some of the 
energy would effectively travel mainly in the 
high velocity strata and give rise to an apparent 
uniform high velocity, parallel to the bedding 
plane. Across the bedding plane, the thin streaks 
of high velocity would have an effect proportional 
only to their total thickness and the average 
would be approximately the low velocity of the 
thicker layers of the medium. 


One may therefore still retain the assumption 
of isotropy with regard to the individual strata 
and attribute the apparent anisotropy to the 
high velocity refracted rays that will arrive first 
at the recorder, even if the source be placed in 
one of the lower velocity strata. 

Returning now to the general theory of refrac- 
tion shooting it may be concluded that the 
“refraction” paths usually assumed in geophys- 
ical analyses are not only geometrically necessary 
from the point of view of the observed time 
distance curves, but are in fact a natural accom- 
paniment of the wave interference phenomena 
arising from the reflection of spherical elastic 
waves from interfaces between media of different 
elastic constants. 

The writer is indebted to Mr. R. D. Wyckoff 
and to Drs. L. J. Peters, E. A. Eckhardt and P. 
D. Foote for helpful discussion and criticism. 
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The static and the impulse breakdown of paraffin oil 
and xylol were investigated by means of the electro-optical 
shutter. To obtain satisfactory operation of the shutter it 
was found necessary to use a double shutter consisting of 
three crossed nicol prisms and two Kerr cells. Because 
of the liquids being only commercially pure and of the 
occluded air and moisture, several.types of sparks were 
observed for a given method of applying the voltage. 
For the early stages of static breakdown, however, the 


ionization was usually more intense at the anode. With 
impulse breakdown in paraffin oil there was observed no 
increase in uniformity of type of breakdown as compared 
to static breakdown. With xylol, on the other hand, a 
fair degree of uniformity in type was observed. This 
uniformity was due to the extremely short time of voltage 
application, thus causing a decrease in the secondary 
effects resulting from impurities and occluded gasses. 


INTRODUCTION 


HE study of electrical breakdown in liquids 

has been an object of increasing interest for 

the past several years. Undoubtedly the extensive 
use of oil for insulating and for quenching sparks 
has fostered a good deal of this interest. In spite 
of the fact that there has been considerable study 
and experimentation along this line there is very 
little known about the mechanism of breakdown 
in liquids. This is probably due to the complicated 
nature of the breakdown and still more to the 
extremely short times in which it takes place. In 
the electro-optical shutter we find an apparatus 
that is capable of operating in the extremely 
short time intervals required. That is, with this 
shutter it is possible to observe the spark in its 
initial stages of formation. By means of such 
observations on sparks in air F. Dunnington has 
been able to determine the initial regions of 
ionization and how they progress.’ From these 
data he has been able to learn something about 
the space charges. With this shutter Dunnington 
was also able to determine the rate of fall of 
voltage across the gap when a spark occurs. 
Observations are now in progress to determine 
the spectrum of the light emitted by the initial 
stages of breakdown. From these data it is hoped 
to determine the kind of ions that are being 
formed. It is seen from this that the electro- 
optical shutter has already proven itself a 


1 Frank G. Dunnington, Phys. Rev. 38, 1506 (1931). 


powerful tool in the study of electrical break- 
down. 

In the present paper there are presented the 
first results that have been obtained with the 
shutter for the breakdown of liquids. In these 
experiments two methods of applying the voltage 
were used. (1) The voltage was raised to its 
maximum value in about one-half second and 
then held there until breakdown occurred. This 
will be called static breakdown. (2) The voltage 
was applied only for 0.3X10-* second. 
breakdown resulting from this second method 
will be called impulse breakdown. Observations 
of both static and impulse breakdown were made 
on paraffin oil and on xylol. As would be expected 
a manner of formation was observed which is 
decidedly different from that of the ‘ionization 
breakdowns”’ observed in air. 


. 
THE ELectro-OpTicAL SHUTTER 


In developing the electro-optical shutter for 
use in the study of electrical breakdown of 
liquids the writer first attempted to use the 
simple type (see Fig. 1) that was employed 
successfully by F. Dunnington in the study of 
sparks in air. This simple shutter was however 


Fic. 1. Diagram of the electrical and optical circuits of the 
single electro-optical shutter. 
29 
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30 HAROLD W. 
found to be unsatisfactory for sparks in liquids. 
That is, the shutter would not close effectively 
enough to shut out the bright light emitted by 
the later stages of the spark. To eliminate this 
difficulty a double shutter was tried which 
proved to operate satisfactorily. (See Fig. 2.) 
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Fic. 2. Diagram of the electrical and optical circuits of the 
double electro-optical shutter. 


In order to make clear the reasons for the 
improvement of the double shutter over the 
single shutter it will first be necessary to point 
out the characteristics of the electro-optical 
shutter. The operating characteristics of the 
simple shutter have been fully discussed in a 
previous paper,’ and so will only be briefly given 
here. A diagram of the shutter is shown in Fig. 
1. It will be noted that there is an optical 
system and an electrical circuit. The optical 
system is composed of two nicol prisms, a Kerr 
cell, a mirror and lenses. The two nicol prisms are 
crossed so that no light from the spark passes 
through them when the Kerr cell is disconnected. 
The Kerr cell is oriented so that the direction of 
its electrical field is at 45 degrees to the plane of 
polarization of the light passed by the first nicol 
prism. With this arrangement the light from the 
spark can pass through the shutter only when a 
voltage is on the Kerr cell. The transmission of 
* the shutter as a function of the voltage on the 
Kerr cell is shown in Fig. 3. In this diagram 100 
percent voltage is the voltage at which the 
shutter has its maximum opening. 

In the electrical circuit R; is a resistance of 
about 10* ohms and serves to limit the rate of the 
building up of voltage on C, and hence the spark 
gap. The distributed resistance between C, and 
the spark gap is 200 ohms which is large enough 
to make the discharge of C; (560 m.m.f.) into the 
spark aperiodic. The distributed resistance be- 
tween the spark gap and Kerr cell is adjusted to 
give the best operating conditions of the shutter. 
This best value is determined experimentally, as 
will be described presently. 


2H. W. Washburn, Phys. Rev. 39, 688 (1932). 
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The operation of the shutter will now be 
briefly described. The voltage on the spark gap is 
brought up in about one-half second. This voltage 
is of course also on the Kerr cell and the shutter is 
therefore open. There is, however, no light 
passing through the shutter until the spark 
begins to form. At the instant the spark begins to 
form the voltage drops across the gap. Because of 
this drop in voltage the Kerr cell begins to 
discharge. When it has discharged to an extent 
which makes the light transmitted by the shutter 
too feeble to be observed, the shutter is effectively 
closed. Consequently the light from the spark 
which is observed is that which passes through 
the Kerr cell up to the time at which the shutter 
effectively closes. The stage of development that 
is observed can be altered either by changing the 
rate of discharge of the Kerr cell by appropriate 
changes in the Kerr cell leads or by changing the 
length of light path between the gap and the 
shutter. This change in light path can be 
conveniently obtained by means of a mirror and 
lens as shown in Figs. 1 and 2. In making 
observations on sparks in liquids this change in 
length of light path was accomplished by two 45 
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Fic. 3. Light transmission curves. A, Single shutter. B, 
Double shutter. 


degree mirrors which sent the light back a 
little to one side of the spark chamber. This was 
done to avoid unnecessary distortion of the 
image which would result from again passing the 
light through the imperfect lens surface formed 
by the plane windows and adjacent liquid in the 
spark chamber. 

The problem of adjusting the Kerr cell spark 
gap circuit for satisfactory operation is not as 
simple as it first appears. It is not only necessary 
that the Kerr cell discharge rapidly enough but 
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also that it remain nearly discharged until the 
spark is quenched. The Kerr cell leads have 
inductance as well as resistance and so the 
discharge of the Kerr cell is of an oscillatory 
nature. If the resistance of the Kerr cell leads is 
too small the oscillations of voltage on the cell 
will be large enough to effectively re-open the 
shutter. If this resistance is too large the Kerr 
cell will discharge too slowly and the closing time 
of the shutter will not be short in comparison 
with the formation time of the spark. Both of 
these effects pass light from the later stages of 
the spark and so it is not a difficult matter 
experimentally to determine the resistance which 
gives the most satisfactory operation of the 
shutter. 

Although it was possible to obtain good results 
with this simple shutter for sparks in air it was 
not found possible to obtain satisfactory oper- 
ation for sparks in paraffin oil or xylol. No 
resistance could be found which would not pass 
so much light from the later stages of the spark 
that the earlier stages were either masked or 
obliterated. This was probably due either to the 


later stages of the sparks in liquids being so’ 


much brighter than earlier stages, or to a 
relatively high voltage appearing across the gap 
after the spark had formed. Using the cathode- 
ray oscillograph, Rogowski has observed that in 
oil the voltage may return to a considerable value 
(about 1 kv for a 0.15 mm gap) in about 10-° 
seconds after the first major drop at the beginning 
of breakdown. 

By employing a double shutter as shown in 
Fig. 2, both of the above-mentioned difficulties 
can be overcome. The double shutter consists 
essentially of two single shutters optically in 
series and electrically in parallel. The second 
nicol prism is crossed with the first and the third 
with the second. Both Kerr cells are oriented so 
that their electrical field is at an angle of 45 
degrees to the plane of polarization of the 
incident light. The transmission curve of this 
double shutter is shown in Fig. 3. It is obtained 
by squaring the ordinates of the transmission 
curve for the single shutter. It will be noted that 
the double shutter has the same transmission at 
100 percent voltage but that its transmission for 
the lower voltages is considerably lower than 


3 W. Rogowski, Arch. f. Elektrotechnik 18, 510 (1927). 


that of the single shutter. As regards the rapidity 
of closing of the double shutter it is seen that its 
transmission falls more rapidly with voltage. On 
the other hand the rate at which the voltage 
drops is reduced since the Kerr cell capacity has 
been doubled. Calculations show that these two 
factors just about neutralize each other. Conse- 
quently the rate of closing of the two shutters is 
approximately the same. The double shutter, 
however, has the disadvantage of having a much | 
longer light path which reduces the amount of 
light that can be collected from the spark and 
passed through to the observer. The double 
shutter would therefore be used only where a 
more complete closing is necessary to prevent the 
passage of light from the later stages of the 
spark. 


THE ELECTRICAL BREAKDOWN OF LIQUIDS 


Results of other experiments*: on the break- 
down of liquid dielectrics have shown that the 
breakdown voltage for a given temperature and 
pressure varies considerably with the duration 
of the application of voltage and with the purity 
and gas content of the liquid. Nikuradse points 
out that, with transformer oil which has been 
carefully purified and outgassed, the breakdown 
voltage is extremely high, being of the order of 
1.310° volts per centimeter. Experiments on 
thin films indicate that the breakdown in pure 
outgassed oil is a pure “ionization” breakdown.’ 
A pure ionization breakdown is conceived as an 
ionization by collision process that results from 
the destruction of the dielectric properties of the 
liquid as distinguished from the destruction of 
these properties by secondary effects such as 
ionizations caused by the presence of impurities, 
occluded air or moisture, or by a temperature 
rise resulting from the conduction current. In 
the article just referred to Nikuradse also shows 
that, as the time of voltage application decreases, 
the breakdown voltage increases and the second- 
ary processes have less effect. For an impulse 
voltage as short as 10~* to 10" second he shows 
that the secondary effects would have completely 


4A. Nikuradse, Arch. f. Elektrotechnik 25, 826 (1931). 

5L. Inge and A. Walther, Arch. f. Elektrotechnik 23, 
279 (1930). 

®H. Edler, Arch. f. Elektrotechnik 25, 447 (1931). 

7A. Joffe, Zeits. Sowj. 1, 255 (1932). 
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vanished and that in all probability a pure 
ionization breakdown would be obtained. 

From the above discussion it appears that one 
should start experimenting with a pure outgassed 
liquid, since the pure ionization breakdown is 
much simpler than the breakdowns resulting 
from many secondary effects. With the simpler 
phenomenon one would have a better chance of 
unravelling the mechanism. On the other hand, a 
single breakdown in the pure outgassed liquid so 
changes its properties that subsequent break- 
downs occur at a much lower voltage. This 
means, of course, that only one observation can 
be made with each sample that is prepared. 
Since in the present experiment the electro- 
optical shutter was being developed for use with 
breakdown in liquids and since this development 
required hundreds of observations, it was not 
practical to start with pure liquids. With the 
knowledge gained from the present experi- 
mentation, however, it is believed that it will not 
be difficult to obtain observations on pure liquids. 
This is to be the next experimental study by the 
writer. 

In the present study observations were taken 
on commercially pure paraffin oil and xylol 
both for static breakdown and an impulse break- 
down of 0.3 X 10~* second. It was necessary to use 
an impulse voltage as long as this since it takes 
between 10-* and 10~’ second to charge up the 
Kerr cells. This will be made clear in the dis- 
cussion of the impulse circuit which will be given 
in the next section of this article. The results 
obtained with paraffin oil and impulse break- 
down were nearly as complex as those obtained 
with static breakdown. On the other hand the 
impulse breakdown of xylol gave much more 
uniform results than were obtained with static 
breakdown. 


Tue IMPULSE CIRCUIT 


The impulse circuit employed makes use of the 
usual scheme of initiating the impulse by means 
of a secondary spark gap. The circuits ordinarily 
employed, however, were not directly applicable 
as they require a very small capacity of the gap 


to which the impulse is applied. In our case the. 


Kerr cells are connected across this gap and a 
small capacity is therefore not realized. An 
additional difficulty is that the spark gap-Kerr 
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cell circuit is oscillatory. It is desirable that the 
impulse voltage be applied in such a way that the 
Kerr cells and spark gap come up to a given 
voltage and then remain there. If oscillations 
occur this condition of course will not be ob- 
tained. It was approximately obtained by placing 
a large resistance between the gap and the supply 
condenser (R, in Fig. 4). 


Fic. 4. Diagram of the impulse circuit used in connection 
with the electro-optical shutter. 


The operation of this impulse circuit is made 
clear by a study of the diagram shown in Fig. 4. 


“The voltage is applied to the transformer which 
brings up the voltage on Gs, C2, and C, in about - 


one second. This rate of voltage build-up is 
controlled by the filament current in Kz and by 
R, which is about two megohms. The voltage is 
raised until a spark occurs at G2. This spark may 
be thought of as suddenly connecting C, to 
the test gap G, and Kerr cells through Ry. The 
resistance R, was 200 ohms which is large enough 
to reduce the oscillations to less than five percent 
of the voltage on C,. A larger resistance would 
make the discharge at G, too feeble. The re- 
sistances R, were 40 ohms apiece which was 
sufficient to make the spark at G2 last somewhat 
longer than 0.3X10-* second. The inductance 
and resistance L; and R; were adjusted so as to 
limit the time to which full voltage was applied 
to G, to about 0.3 10~ second. The calculated 
test gap voltage-time characteristic is the solid 
curve shown in Fig. 5. The dotted curve is the 
breakdown voltage of a 0.8 mm gap in unpurified 
transformer oil plotted as a function of the time 
of voltage application. This curve was taken 
from curves given by R. Naher.® It is seen that 
the voltage-time characteristic of oil is such that 


8 R. Naher, Arch. f. Elektrotechnik 21, 169 (1929). 
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Fic. 5. Curve A, calculated impulse voltage applied to 
test gap. Curve B, impulse breakdown voltage character- 
istic of unpurified transformer oil. 


an impulse voltage with a long trailing foot, such 
as the one shown, can be used without danger of 
incurring any appreciable variation in the time or 
the voltage of breakdown. 


EXPERIMENTAL PROCEDURE 


In taking observations of static breakdown the 
potential was first adjusted to 45,000 volts. This 
was the voltage that gave maximum opening of 
the shutter. The voltmeter used consisted of 
fifty feet of one-quarter inch glass tubing 
through which distilled water was circulated. 
The current taken by this water column was 
about 100 microamperes at 50,000 volts. There 
was therefore no appreciable ripple introduced by 
the voltmeter. The resistance of the water column 
was determined each day with a 5000 volt 
(calibrated) electrostatic voltmeter. After the 
voltage was set, the gap was adjusted until the 
breakdown occurred in one to ten seconds after 
the application of voltage. About ten seconds 
were allowed between each breakdown in order 
that the gas bubbles formed by the spark could 
rise out of the field between the electrodes. As 
was previously stated, liquids of commercial 
purity were used. With this purity about fifty 
breakdowns could be observed on the same 
sample before any noticeable effect on the break- 
down voltage occurred. The electrodes, which 
were eight millimeter brass spheres, were cleaned 
and polished after every 200 breakdowns. With 
the small energies of discharge employed there 
was very little corrosion of the electrodes. 

In taking observations of impulse breakdown, 
the spark gap Gz was adjusted so that 35,000 


volts would occur across the test gap G:. The 


test gap electrodes were then set at the maximum 
spacing at which a spark would occur. The rest 
of the procedure was similar to that just given 
for static breakdown. 


DETERMINATION OF THE OPEN TIME OF 
THE SHUTTER 


The length of time that the shutter remains 
open after the beginning of breakdown can be 
calculated from the circuit constants and a 
knowledge of the rate at which the gap voltage 
falls.* For our present qualitative study, however, 
an accurate knowledge of this time is not 
necessary. Indeed it would be impossible to 
make an accurate calculation without an accurate 
knowledge of the rates of voltage drop for our 
particular conditions. A quantity known as the 
“‘path time’”’ will therefore be used as a measure 
of the open time of the shutter. The path time 
is defined by the following expression: Path time 
=T,—T i, where T,,=length of wire path be- 
tween the spark gap and Kerr cell divided by the 
velocity of light; 7,=length of light path be- 
tween the spark gap and Kerr cell divided by the 
velocity of light. An inspection of Fig. 2 will 
show that this path time is the interval of time 
between the instant at which the light from the 
beginning of the spark reaches the Kerr cell and 
the instant at which the beginning of the voltage 
wave from the spark gap reaches the Kerr cell. 
The actual open time of this shutter is somewhat 
longer than the path time.': ? If the time that it 
takes the voltage to drop for our conditions is 
10-* second, about 0.810-* second should be 
added to the path time to obtain the open time of 
the shutter. This value was obtained by making 
some calculations for the particular leads used 
and by assuming the voltage to drop in 10% 
second. 


EXPERIMENTAL RESULTS 


a. Static breakdown of paraffin oil 


The observations of the static breakdown of 
paraffin oil were taken with a potential of 45,000 
volts and the gap adjusted for breakdown to 
occur in from 1 to 10 seconds after the application 
of voltage. This gap was found to be 2.1 milli- 
meters. Observations of the earliest part of the 
breakdowns were made visually as the light 
emitted was too feeble to photograph. With the 
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shutter set for a path time of —0.8 X 10-* second 
no light at all could be seen. This is to be expected 
since an approximate correction of this path time 
by the figure given in the preceding section gives 
zero time. In other words, with this setting, the 
shutter closes at the instant the spark begins to 
form. For a setting of —0.610-% second path 
time the spark is visible. Due to the commercial 
purity of the oil several types of sparks were 
observed. They all have, however, certain com- 
mon characteristics, as follows. As soon as the 
spark is visible at all it is visible throughout the 
greater part of the distance between the elec- 
trodes. Some parts are brighter than others but 
completely dark spots usually occur singly or in 
pairs and only cover a relatively small portion of 
the gap. Occasionally an exception to this rule 
occurs when the spark begins by emitting light 
only at several points in a line between the 
electrodes (No. 7, Fig. 6). Fig. 6 shows a few 
freehand drawings which were made of each 
spark at the time of its observation. 


Fic. 6. Freehand drawings of initial stage of static break- 
down of paraffin oil. Path time = —0.6 x 10-8 second. 


The manner in which the spark in oil begins 
differs markedly from that of sparks in air. As 
was found by Dunnington the sparks in air 
start by intense ionizations at the electrodes and 
under certain conditions with an additional 
region between the electrodes. 

For path times of 0.06 10-%, 0.3 10-5, and 
0.8 second the light intensity was suffi- 
ciently bright to be photographed. A few of the 
pictures obtained are shown in Fig. 7. The 
numbers under the sparks indicate the number of 
times that type of spark occurred out of fifty 
sparks taken with the same shutter setting. It 
will be noted that for the shortest path times 
photographed (0.06X10-* sec.), the sparks 
having one or two points of intense ionization at 
the anode occur most frequently (types 1 and 2 in 
Fig. 7). At the two later times (0.310-* and 
0.8 X 10~* sec.) photographed, the predominating 


®F, Dunnington, Phys. Rev. 38, 1535 (1931). 
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sparks are those having either several points ot 
intense ionization or nearly uniform ionization 
throughout their length. These types are shown 
in columns 9 to 14 in Fig. 7. 

The pictures in the row labeled ‘“‘no shutter” 
were taken with the nicol prisms uncrossed and 
the Kerr cells disconnected. The diaphram in the 
first lens was cut down so that the brightest parts 
would give densities well below the saturation 
point of the plate. These pictures therefore show 
an integration of the intensities emitted through- 
out the duration of the entire spark. The length 
of this time is of the order of magnitude of 10~° 
second. It is seen that the regions of intense 
ionization in over half the cases are differently 
located than in the case of any observed spark at 
the early stages of breakdown. The most pre- 
dominate type apparently has its highest in- 
tensity at both electrodes. It would be well to 
mention at this point that the intensities of light 
emitted by the spark are decidedly not pro- 
portional to the density of blackness on the 
pictures shown. The contrast is exaggerated by 
the photographic process. This explains why 
some points occur in the spark, even with no 
shutter action, that are perfectly white. 


b. Impulse breakdown of paraffin oil 


In taking observations of the impulse break- 
down of paraffin oil the voltage was adjusted to 
35,000 and the gap set at the maximum spacing 
at which a spark would occur. The spacing of the 
electrodes thus obtained was about 0.4 milli- 
meter. The observations of all the impulse 
breakdowns were made visually as the images 
were too weak to photograph. No light from the 
spark was visible with a path time as short as 
—0.8X10-* second as was the case with static 
breakdown. With the shutter set for a path time 
of —0.4X10-* second a very faint haze is seen 
between the electrodes. This haze is not uniform 
but takes various shapes. With a path time of 
—0.3X10- second bright spots begin to appear 
in this haze. These bright spots apparently do 
not appear in any particular region but occur as 
often somewhere between the electrodes as at 
the electrodes. For path times as long as 0.5 X 10~° 
second about eight out of ten sparks show nearly 
uniform ionization throughout their length. Fig. 
8 shows some freehand drawings made of these 
sparks. 
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TYPE NO. 


PATH TIME | Tet iL 
0.06x10-° SEC. |-| | | 
1|2le 
PATH TIME | 7 | ta 
0.3x10°° SEC. | —| _ 
de 
PATH TIME |} | 7 
— — | |.—._ | - 
0.8x10°8 SEC. | 
9 si2ii4 8 |4 18 


Fic. 7. Photographs of single sparks obtained by the static breakdown of paraffin oil. 


It will be noted that, even though the break- 
down voltage to impulse potentials for this oil is 
not far below that of pure outgassed oil, the 
manner of breaking down is not uniform as would 
be expected of pure outgassed oil. 


c. Static breakdown of xylol 


The observations on xylol were made for a 
potential of 45,000 volts and the gap, adjusted as 
before, was found to be 2.3 millimeters. With the 
apparatus adjusted for a path time of —0.8 x 107° 
second nothing could be observed. This agrees 
with the results obtained with paraffin oil. At a 
path time of —0.6X 10-* second only a very dim 
haze between the electrodes was visible. At a path 
time of —0.410-> second regions of intense 
ionization began to appear, usually at the anode. 


Sometimes the spark would be fairly bright 
throughout its length. As these observations were 
all visual it was necessary to record the results by 
freehand drawings. The types observed at this 
early time were, however, the same as those 
recorded photographically at a path time of 
0.06 X 10-* second (see Fig. 9). The drawings are 
therefore not shown. The essential difference 
between the beginning of the spark in xylol and 
in paraffin oil is that the ionization in xylol at 
first is more localized and larger gaps in the 
spark appear. 

For path times of 0.0610-°, and 
0.8 second the light intensity was bright 
enough to photograph. The several types of 
sparks obtained are shown in Fig. 9. As before, 
the numbers under the pictures indicate the 
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Fic. 8. Freehand drawings of sparks obtained by the impulse breakdown of paraffin oil. 
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TYPE NO. 8] 9 

+ 

PATH TIME 
0.06x10°° SEC. | — 

5 3/1 

SEC. | 

a8] 3] 2/2 

0.8x10-°* SEC, | 

16} 9/10 

NO SHUTTER | — 

6 8 is|}6]|3 8 4 


Fic. 9. Photographs of single sparks obtained by the static breakdown of xylol. 


number of times that spark occurred out of a 
group of fifty sparks. It will be noted that for the 
three path times shown the type of spark oc- 
curring most frequently has a single or double 
region of intense ionization at the anode. It will . 
be recalled that this was the case for the earliest 
stage photographed with paraffin oil. 

The column labeled ‘‘no shutter’’ was obtained 
in the same way as for paraffin oil. It will be 
noted that sparks having a double region of 
intense ionization at the cathode occur twice as 
often as those having this region at the anode. 
Another great difference between these sparks 
and those obtained for the initial stages of 
breakdown is that sparks having regions of 
intense ionization at both electrodes are quite 
numerous whereas for the early stages this type 
very seldom occurs. 


d. Impulse breakdown of xylol ‘ 


The impulse breakdown of xylol was observed 
for an applied potential of 35,000 volts. The gap 
was adjusted as previously described and found 
to be 0.5 millimeter. Again with the apparatus 
set for a path time of —0.8 X 10~ second no light 
from the spark could be observed. At a path 
time of —0.3 x 10-* second the spark appeared as 
a faint haze between the electrodes and occasion- 
ally in addition a bright spot at just one or at 


both electrodes (see Fig. 10). At a path time of 
zero seconds the spark was usually fairly bright 
all the way across with the regions at one or both 
electrodes brighter, as is shown in Fig. 10. A 
bright region, however, occurred at the anode 
more often than it did at the cathode. At a path 
time of 0.510-* second the spark was always 
bright throughout its length. Sometimes a region 
brighter than the rest would occur at one or both 
electrodes but more often at the anode. It is 
interesting to note here that, with impulse 
breakdown of commercially pure xylol, fairly 
uniform results were obtained. 


DISCUSSION 


From the experimental results just described 
it is seen that there are several types of sparks 
which occur in commercially pure paraffin oil and 
xylol. In fact the results are so complex that it is 
difficult to draw any definite conclusions as to the 
mechanism of breakdown. It might be instructive, 
however, to suggest an explanation for the 
intense ionization observed at the anode in the 
early stages of breakdown. From a consideration 
of the views on cataphoresis of Lenard and of 
Coehn and from the experiments of Dyk," it is 
the belief of L. B. Loeb that when liquids, such as 


© Karl Dyk, Phys. Rev. 31, 913 (1928)-: 
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those studied by the writer, are subjected to a 
high electric field the gas bubbles present in the 
liquid will move to the anode. Since in the present 
experimentation no attempt was made to remove 
the gas from the liquids, an accumulation of gas 
bubbles at the anode would be expected soon 
after each application of voltage. These gas 
bubbles might be expanded by the heating and 
consequent evaporation of the surrounding liquid 
into the bubbles as suggested by Inge and 
Walther.’ With a sufficient expansion the break- 
down would start by ionization within the 
bubble. On the other hand the mere accumulation 
of the gas bubbles at the anode might be suffi- 
cient to initiate the breakdown. In either case 


+ + + + + + 
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Fic. 10. Freehand drawings of single sparks obtained by the 
impulse breakdown of xylol. 


their presence at the anode would in part 
account for the intense ionization observed in 
this region during the early stages of breakdown. 
It is of interest to note further that the time 
required to move the bubbles to the electrode 
would explain the increase in breakdown voltage 
with the increase in frequency of the applied 
voltage (in the range zero to 800 cycles) as 
observed by Draeger."' 

From the results of impulse breakdown the 
following observation is of interest. The fact that 
uniform results were obtained for the impulse 
breakdown of xylol indicates that the secondary 
effects have been nearly eliminated by the short 
time of application of voltage. On the other hand 
the results for paraffin oil indicate that the 
secondary effects in its breakdown have not been 
eliminated. Inge and Walther came to this same 


" K. Draeger, Arch. f. Elektrotechnik 13, 366 (1924). 


conclusion from their experiments. They found 
that the breakdown voltage of xylol for an 
impulse of 0.210-* second did not vary with 
pressure, temperature, or occluded air and varied 
only slightly with impurities. On the other hand 
the impulse breakdown of transformer oil showed 
a variation with pressure and temperature. As 
was previously stated Nikuradset showed that 
these secondary effects could be eliminated in the 
breakdown of transformer oil for impulses as 
short as to second. 

It should be emphasized that it is very de- 
sirable to obtain pictures of sparks in very pure 
outgassed liquids. When this has been done the 
effect of impurities and occluded air can be 
experimentally determined. Such experiments 
would throw light upon the present results, and a 
better understanding of the mechanism of break- 
down would undoubtedly result. It is hoped that 
it will be possible in the near future to obtain 
results along this line. 


CONCLUSIONS 


(1) In the static breakdown of paraffin oil 
ionization occurs throughout almost the entire 
distance between the electrodes in 2X10~° 
second. This ionization is not uniform but may 
have one, two, or sometimes several regions of 
higher intensity. At the earliest time at which the 
spark can be photographed it was definitely 
determined that these regions of high intensity 
usually occur at the anode. 

(2) In the static breakdown of xylol the 
ionization in the initial stages of breakdown is 
more localized than with paraffin oil. At the 
earliest time that the spark could be clearly seen 
(T,= —0.4X10- sec.) the type most often ob- 
served had one or two regions of intense ioniza- 
tion at the anode. 

(3) By applying an impulse voltage of 
0.3X10-* second to paraffin oil no gain in 
uniformity of the type of spark was obtained. 

(4) The results for the impulse breakdown of 
xylol showed an improvement in uniformity over 
the static breakdown. The early stages of the 
spark showed an intense ionization at one or both 
electrodes. 

In conclusion I wish to express my appreciation 
to Professor E. O. Lawrence for his valuable 
suggestions and cooperation in carrying out this 
experimentation. 
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